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ABSTRACT 

Local heat transfer and pressure coefficients around 
a right circular cylinder were experimentally determined at 
Reynolds numbers ranging from 57,000 to ^95,000. The turbu- 
lence intensity of the free stream was approximately 0.5%. 
The cylinder was heated externally with a constant heat flux 
and simultaneously cooled in a cross flow of air. A uniform 
heat source was provided by energizing Nichrome ribbon with 
a constant, measured current. The steady state surface tem- 
perature increase above the air temperature was indicated by 
thermocouple information. Employing Newton’s law of cool- 
ing, heat transfer coefficients were determined as a func- 
tion of the angular location from the forward stagnation 
line. The free stream dynamic pressure and surface pressure 
distribution around the cylinder were obtained using static 
pressure pickups and their associated manometers. A com- 
parison is made between the heat transfer and pressure data 
collected. Experimental results compared within about six 
per cent of the theoretical solution of Schuh. 
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NOMENCLATURE 



English Letter Symbols 



Ai Cross-sectional area of wind tunnel 

in settling chamber 

A 2 Cross-sectional area of wind tunnel 

at entrance to test section 

A Cross-sectional, area of Nichrome 

n ribbon 



A s 

C 

P 



D 

e 



f 

F 

h 

h e 



i 

k 



Nu 

Pi 



External surface area of an ele- 
ment of Nichrome ribbon 

Pressure coefficient as a function 
of angular location on cylinder 
surface 

Diameter of cylinder 

Root mean square voltage of hot 
wire anemometer 

subscript; Quantity evaluated at 
film temperature 

Proessling number, F = Nu/v / Re 

Convective heat transfer 
coefficient 

Local convective heat transfer 
coefficient as a function of 
angular location 

Current 

Thermal conductivity of air 
evaluated at film 
temperature 

Thermal conductivity of Nichrome 

Nusselt number, Nu = hD/k 

Static pressure of air in 
settling chamber 



ft 2 

ft 2 

ft 2 

ft 2 



ft 

volts 



BTU 

hr-ft 2 -°F 
BTU 

hr-ft 2 -°F 

amps 

BTU 

hr-f t-°F 
BTU 

hr-ft-°F 

lbf/ft 2 
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English Letter Symbols 



P 2 


Static pressure of air at 
entrance to test section 


lb f /ft 2 


P 3 


Static pressure of air in test 
section 


lbf/ft 2 


p e 


Surface pressure on cylinder as 
a function of angular location 


in of water 


Pr 


Prandtl number 




APl-2 


Static pressure drop between 
settling chamber and entrance 
to test section 


cm of water 


Ap 

cylinder 


Pressure difference between sur- 
face pressure on cylinder and ' 
static pressure at entrance to 

test section, Ap = P Q -P 2 

* cylinder 0 


cm of water 


q 


Heat flux 


BTU/hr-ft 2 


q e 


Local heat flux from an element 
at some angular location 


BTU/hr-ft 2 


q i 


Dynamic pressure at any location, i 


lbf/ft 2 


q ideal 


Ideal dynamic pressure in test 
section 


lbf/ft 2 


q actual 


Actual dynamic pressure in test 
section 


lbf/ft 2 


Q 


Heat transfer rate 


BTU/hr 


Q k - 


Heat transfer rate through 

Nichrome element by conduction 


BTU/hr 




Heat generation rate within any 
Nichrome element, i 


BTU/hr 


Q h 


Heat transfer rate by convection 
from external surface of 
Nichrome element 


BTU/hr 


S?ad 


Heat transfer rate by radiation 
from external surface of 
Nichrome element 


BTU/hr 


Re 


Reynolds number. Re = UD/v 


• 
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English Letter Symbols 



R 

T 



i 

AT 

Tu 

U 

Ux 

U 2 

u. 



u 

/ u* 1 " 
V - 

V 0 

Ax 

X 



Resistance of main heating circuit 

Surface temperature of cylinder 

Local surface temperature of 
cylinder at some angular 
location, 0 

Temperature of ambient air in wind 
tunnel 

Surface Temperature at any element 
of Nichrome ribbon, i 

Local temperature difference 

between cylinder surface and the 
ambient air, AT = T -T 

0 oo 

Turbulence intensity level 

Free stream velocity 

Free stream velocity in settling 
chamber 

Free stream velocity at entrance 
to test section 

Free stream velocity in test 
section 

Mean free stream velocity 

Mean fluctuation in free stream 
velocity 

Voltage in heating circuits 

Initial probe voltage of hot wire 
anemometer at zero free stream 
velocity 

Arc length of an element of 
Nichrome ribbon 

Tunnel Correction Factor 



ohms 

op 

°F 

°F 

°F 

op 

ft/sec 

ft/sec 

ft/sec 

ft/sec 

ft/sec 
ft /sec 

volts 

volts 

ft 
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Greek Letter Symbols 



sb 



'wb 



Pi 

P 2 



Thermal emissivity of Nichrome 

Wind tunnel boundary correction 
for solid blockage 

Wind tunnel boundary correction 
for wake blockage 

Total wind tunnel boundary correc- 
tlon > e t = e sb +e wb 

Density of air in settling chamber 

Density of air at entrance to test 
section 



Pa 



Density of air in test section 



o)^ Uncertainty interval of quantity i 

v Kinematic viscosity of air evaluated 

at film temperature 



lbm/ft 3 
lbm/ft 3 

lbm/ft 3 
ft 2 /hr 
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I. INTRODUCTION 



Local heat transfer and pressure coefficients around a 
cylinder were experimentally determined at ten different 
Reynolds numbers that ranged from 57,000 to 495,000. The 
turbulence intensity level of the free stream was determined 
with a hot wire anemometer and found to be approximately 
0.5 per cent. 

The cylinder was heated externally with a constant heat 
flux, and simultaneously cooled in a crossflow of air by . 
forced convection heat transfer. 

A uniform surface heat flux was provided by helically 
winding Nichrome ribbon into shallow surface grooves that 
extended around the circumference of the cylinder. Approx- 
imately half of the thirty-two inch length of the cylinder 
was wrapped with the main heating circuit. Two smaller 
guard heater circuits were similarly installed on each end 
of the main circuit. Thermal energy dissipated by Joulean 
heating, as described by Wood in Ref. 1, was then applied to 
the surface by energizing the electrically resistant ribbon 
with a constant, measured current. The guard heater circuits 
were separately controlled. 

The 4.5 inch diameter cylinder spanned the entire height 
of the wind tunnel test section, permitting a two dimension- 
al study. To prevent any internal free convection, the 
cylinder was filled with foam rubber insulation. 
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After the thermal capacity of the cylinder was saturat- 
ed, a steady state condition was achieved. Then all of the 
heat added was removed by convection, less any losses due 
to thermal radiation to the tunnel walls or circumferential 
conduction through the Nichrome ribbon. The guard heaters 
barred heat from being conducted out the ends. The local 
surface temperature increase above the ambient air tempera- 
ture, delta T, was indicated by copper-constantan thermo- 
couples welded to the inner surface of the Nichrome ribbon. 
The test section, which was located at the axial midpoint of 
the cylinder, was instrumented with four of these thermo- 
couples welded to the center wrap of the Nichrome ribbon at 
ninety degree intervals. Since the cylinder was mounted on 
a turntable, delta T could be measured at any desired angular 
position. Employing Newton’s law of cooling, local heat 
transfer coefficients were then determined as a function of 
angular location from the forward stagnation line. 

Two static pressure taps were also installed on the 
cylinder, ninety degrees apart. These were connected to 
U-tube manometers which were referenced against the static 
pressure. Local surface pressure variations around the 
circumference of the cylinder could then be measured direct- 
ly. The free stream dynamic pressure was measured by two 
static pressure pickups installed in the wind tunnel, and 
indicated on the "air speed" manometer. The measurement 
system, described in detail in Section III, provided 
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sufficient information to evaluate the pressure coefficient 
at any desired angular position. 

Local experimentally determined heat transfer coeffi- 
cients were compared with the measured pressure coeffi- 
cients in order to correlate heat transfer trends with the 
various flow regimes which developed around the circumfer- 
ence of the cylinder. The location of the laminar, transi- 
tional and turbulent boundary layers as well as separation 
and its attendant wake were found, as expected, to produce 
definite trends in the heat transfer data. 

The Froessling number, Nu/i/Re", was determined at various 
angular locations and was compared with the experimental 
and theoretical results of other researchers. The low free 
stream turbulence intensity maintained throughout the pres- 
ent investigation led to experimental results that were 
significantly lower than those previously reported. Experi- 
mental results were found to agree within about six per 
cent of the theoretical model developed by Schuh for predict- 
ing local heat transfer coefficients on the forward portion 
of a -constant heat flux cylinder. 



II. BACKGROUND 



Almost immediately after beginning a study of heat 
transfer from a cylinder, two absolute facts become appar- 
ent. First comes a total realization that the heat trans- 
fer mechanism is closely interrelated with the associated 
fluid mechanics. An adequate treatment of the former, 
whether it be analytical or experimental in nature, requires 
a basic understanding of the movement, pressure, and general 
behaviour of the fluid particles which conduct and convect 
thermal energy. 

Secondly, it is certainly true that a vast amount of 
research has been devoted to each aspect of the thermal- 
hydrodynamic interface around a heated cylinder placed in a 
uniform stream. • A mere listing of all that has been pub- 
lished would be voluminous. Instead the approach has been 
taken to review the better known, standard works and some 
that are unique: The details of the many factors involved 

are discussed within the framework of this review. 

A. FLUID FLOW AROUND A HEATED CYLINDER 

A brief description of this flow model is summarized 
in this section based on discussions contained in Refs. 2-5. 

The potential flow of an inviscid, ideal fluid around 
the surface of a cylinder consists of an acceleration from 
zero velocity at a forward stagnation point, to a maximum 
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velocity at ninety degrees on both sides, and finally a 
deceleration to zero velocity at the rear stagnation point. 
The corresponding surface pressure distribution around the 
cylinder for ideal flow is simply an inversion of the veloc- 
ity distribution. The pressure maxima are located at the 
stagnation points and the pressure minima are located at 
ninety degrees. 

The flow of a real fluid around a cylinder is best 
analyzed in terms of changes that occur within the boundary 
layer that forms near the surface. This thin film of fluid 
which grows from a minimum thickness near the forward stag- 
nation line is a result of the adherence of a viscous fluid 
to the surface, where a no-slip condition of zero velocity 
exists. As long as this layer is attached to a heated sur- 
face it is the medium into which thermal energy is first 
conducted. The fluid motion in the boundary layer then 
convects this heated fluid away from the surface. The 
velocity profile through the boundary layer varies as a 
function of the angular location from the forward stagnation 
line. Hence, the heat transfer rate from the cylinder sur- 
face becomes a function of the angular location. As this 
angular position increases, the boundary layer grows thicker, 
the fluid velocity is increasingly retarded and less heat is 
removed from the surface. 

A boundary layer may be classified as either laminar 
or turbulent . Flow within a laminar boundary layer is char- 
acterized by fluid particles moving in a smooth, orderly 
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fashion and in layers. A turbulent boundary layer is char- 
acterized by individual fluid "packets" flowing in irregular 
patterns within the bounds of a composite motion in the 
direction of flow. A, turbulent boundary layer is conducive 
to a higher heat transfer rate. 

As the kinetic energy of the free stream attenuates 
within the slower moving boundary layer, the pressure near 
the surface becomes less favorable to the established flow 
direction. Finally, when the remaining kinetic energy of 
the fluid near the surface cannot overcome the adverse pres- 
sure rise, the boundary layer separates, as a reverse flow 
is initiated near the surface. Eddies of turbulence result 
wherein packets of fluid "scrub" the heated surface, signif- 
icantly increasing the heat transfer rate. The resultant 
wake behind a cylinder is characterized by low pressure 
relative to the pressure on the forward half of a cylinder. 
At Reynolds numbers ranging from about 60 to 5000 [5], a 
regular pattern of shed vortices persists in the wake; it 
is known as the Karman vortex-street, named after the 
scientist Theodore von Karman. At the higher Reynolds num- 
bers pertaining to this study the wake is characterized by 
irregular turbulent flow. 

Form or pressure drag is a direct result of this low 
pressure region in the wake -caused by boundary layer separa- 
tion. For the range of flow velocities considered in the 
present investigation (50 ,000<Re<500 ,000) the form drag is 
large in comparison to skin friction drag. At higher 
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velocities transition to turbulence occurs prior to separa- 
tion. The separation point moves farther back around the 
cylinder and form drag decreases. 

The characteristics of the boundary layer, the location 
of the separation point, and the characteristics of the 
wake and resultant pressure distribution are dependent upon 
the Reynolds number, 

Rg y > C 1 ) 

where U is the free stream velocity, D is the diameter of 
the cylinder, and v is the kinematic viscosity of the fluid. 
The location of the separation point also depends upon the 
magnitude of the free stream turbulence and surface 
roughness . 

Instability first occurs at the rear stagnation point 
at a Reynolds number of about ten, and half of the total 
drag is then form drag. 

The Karman vortex-street first occurs at Reynolds num- 
bers of about 100. Form drag is then predominant. 

At Reynolds numbers between 1,000 and 100,000 a laminar 
boundary layer separates from the surface at an angular lo- 
cation of 80-85 deg from the forward stagnation line. 

Skin friction drag is negligible in this range. 

Above a Reynolds number of 100,000 the fluid in the 
boundary layer may contain sufficient kinetic energy to 
remain attached for a greater distance and not separate until 
reaching an angle of 130 deg with the flow direction. This 
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was shown experimentally by Fage and Falkner [Refs. 6 and 7] 
to occur, if the boundary layer had become turbulent. The 
critical point of transition from a laminar to a turbulent 
boundary layer moves farther back from the front of a cylin- 
der as the Reynolds number is increased. It generally lies 
in the vicinity of 80-95 deg. Boundary layer transition 
prior to separation is given as a cause of the rapid de- 
crease of the drag coefficient at a Reynolds number between 
200,000 and 500,000. 

At Reynolds numbers above 50,000 the heat transfer down- 
stream of separation is larger than on the front half. This 
difference continues to increase with increasing Reynolds 
numbers . 

Schlichting [5] reviews some previous experimental work 
as follows: Hiemenz measured pressure distribution around a 
cylinder in 1911, determining a separation angle of 81 deg. 
His boundary layer calculations provided a comparable separa- 
tion angle of 82 deg. Extensive experimental data on the 
pressure distribution were later collected by Flachsbart. 

His data exhibited a strong dependency on Reynolds number. 

For subcritical Reynolds numbers he recorded a pressure 
minimum at an angle of about 70 deg and a fairly constant 
pressure over the downstream portion. At Reynolds numbers 
above the critical value he reported a pressure minimum at 
90 deg and generally a closer agreement with potential flow 
theory. Thom investigated the laminar boundary layer on a 
circular cylinder in 1928 at a Reynolds number of 28,000. 
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Page's study [6] of the boundary layer In the region of 
separation from a cylinder included Reynolds numbers from 
100,000 to 350,000. Fage concluded that a critical point 
on a cylinder exists where the boundary layer begins a 
transition to a turbulent state, and that separation occurs 
just beyond the critical point. The critical point was 
shown to be the inflection point in the normal pressure 
distribution curve that immediately follows the point of 
maximum negative pressure. The latter point was also shown 
to move farther around the cylinder with higher Reynolds 
numbers . 

In Ref. 8, Roshko combined free-streamline theory with 
Karman ' s theory of the vortex-street, and developed a solu- 
tion for drag that is dependent on one experimental input, 
the shedding frequency. In the process, he defined a uni- 
versal wake Strouhal number which is a function only of wake 
Reynolds number. 

Roshko 's study of flow past a cylinder at very high 
Reynolds numbers is described in Ref. 9* He observed that 
a second transition" range occurs as the drag coefficient 
"increases in the range, 10 6 <Re<3. 5 X 10 6 , from a value of 
about 0.3 to 0.7, and then levels off at the latter value." 
He also observed that vortex shedding occurs at Re>3*5 x 10 6 . 

Achenbach [10] describes measurements of local pressure 
and skin friction around a cylinder over a range of Reynolds 
numbers from 60,000 to 5,000,000. He determined the loca- 
tions of separation points and transition regions 
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(separation bubbles), and analyzed the data in relation to 
the flow regions involved. 

In Ref. 11, Coder organized and compiled previous ex- 
perimental and theoretical determinations of the boundary 
layer separation point. He also compared this information 
with his own data collected using hot films. His summary 
included separation points at all Reynolds numbers from 
"creeping" flow to "transcritical" flow. His measurements 
were made in the "transition" (or "supercritical") flow 
regime . 

B. THE EVALUATION OF HEAT TRANSFER FROM A CYLINDER 

In the past twenty years it appears that the emphasis 
has shifted from experimental to analytical studies. 
Reference 12 leads one to conclude that about three approx- 
imate solution methods for heat transfer in the forward 
half prior to boundary layer separation have been proposed 
every two years. This trend has evolved a crucial need for 
accurate experimental data to verify the theories. It is 
this motivation which led the author to pursue and refine 
an experimental technique . 

A study of forced convection heat transfer from a cylin- 
der may employ one of two separate approaches. These are 
related to either a constant wall temperature cylinder or a 
constant heat flux cylinder. In the former model, delta T 
Is held constant and a measured heat flux is increased or 
decreased as necessary to maintain a constant temperature 
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wall. In the latter model, the heat flux is held constant 
and the resulting delta T is measured as a function of a 
angular location. 

Both approaches culminate with an application of 
Newton's law of cooling 



q 6 = 



Q 



A 



- V AT) e 



( 2 ) 



where for a specified surface hQ is the local heat trans- 
fer coefficient at angular location 0, Q is the heat 
transfer rate, A is the surface area of interest, q fl is 
the local heat flux and (AT)^ is the local surface tem- 
perature increase above the ambient air temperature. 
Generally, h Q is presented in the dimensionless form of 
the local Nusselt number 



Nu 



V 

k 



( 3 ) 



where D is the diameter of the cylinder and k is the ther- 
mal conductivity of the fluid. 

Experimental and analytical results are often presented 
in the form of the Froessling number, named after N. 
Froessling, which is defined by 



Nu 

/Re 



F 



CM 



1 . Some Previous Experimental Work 



a. Ernst Schmidt and Karl Wenner 

Schmidt and Wenner present in Ref. 13 their 
well known results for heat transfer over the circumference 
of a heated cylinder. A constant surface temperature cyl- 
inder was used in their experiment. Specifically, a hollow 
brass cylinder was maintained at one hundred degrees Centi- 
grade by passing steam vapor through it. A londitudinal 
strip was insulated from the vapor. In fact, the strip 
was physically removed and replaced with a copper element 
that contained an internal electrical heating element. By 
controlling the electrical power into the heating element 
and making use of comparative thermocouple information, the 
strip was maintained at the same temperature as the cylin- 
der wall. Hence, the heat flux was varied to maintain a 
constant wall temperature. A pressure pickup was also in- 
stalled on the cylinder wall. By rotating the entire cyl- 
inder "about a graduated disk," both heat transfer and 
pressure data were collected at various angles with the air 
flow. It was stated that their data were generally only 
taken on half of the cylinder, but symmetry was checked 
with no "substantial difference" between both halves. 

Schmidt and Wenner published the following re- 
sults over a range of Reynolds numbers from 5000 to 426,000: 
"At Reynolds numbers up to around 100,000 high maxima of 
the heat transfer occurred in the forward stagnation point 
at 0 deg and on the rear side at 180 deg, while at around 
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80 deg the heat-transfer coefficient on both sides of the 
cylinder behind the forward stagnation point manifested 
distinct minimums. Two other maximums occurred at around 
115 deg behind the forward stagnation point between 170,000 
and 426,000. At 426,000, the heat transfer at the location 
of these maximums was almost twice as great as in the for- 
ward stagnation point, and the rear half of the cylinder 
diffused about 60 % of the entire heat.... On the rear of 
the cylinder at <}> = 180 deg the heat-transfer coefficients 
for Reynolds numbers of from 25,000 to 426,000 are greater 
than in the forward stagnation point." Heat transfer a- 
round a cylinder was stated to be a strong function of the 
Reynolds number and the angle from the forward stagnation 
line. Additionally, their results for heat transfer at the 
forward stagnation point were said to agree with Squire's 
solution, described in Ref. 14. 

Although Schmidt and Wenner did not directly 
refer to a turbulence factor, Giedt, in Ref. 15, stated that 
there was a minimum of turbulence and that their data were 
taken in a "relatively narrow jet of air with a nearly 
uniform velocity distribution." 
b. R.M. Fand 

In Ref. 16, Fand describes his work which 
involved a constant wall temperature cylinder immersed in 
a uniform stream of water. The Reynolds numbers studied 
were between 11,000 and 64,000. His cylinder was composed 
of a copper core with a surface covering of chromium. 



Thermocouples and electrical heaters were encapsulated in 
hypodermic containers and embedded longitudinally into 
drilled holes in the test cylinder. The high thermal con- 
ductivity of the copper core enabled a uniform surface tem- 
perature to be maintained when the heaters were energized. 

A major design limitation was identified as the low power 
capability of the heaters which provided a temperature in- 
crease of only 4 to 10 °P. Guard heater sections similar 
in design to the test cylinder, were controlled separately 
to prevent any heat flow out the ends. 

Fand’s results were shown to be in close agree- 
ment with the McAdams correlation, which had been based 
upon Reynolds numbers from 0.1 to 200. Fand then concluded 
that this correlation which is used to calculate the aver- 
age Nusselt number, 

(Nu) f = [0.35 + 0 . 56 ( Re ) °* K ] ( Pr )'* (5) 

could be extended to a Reynolds number of 100,000. Pand 
built a strong ' argument for a second preferred form of cor- 
relation formula originated by Douglas and Churchill in 
Ref. 17 , 

0.50 

(Nu) f = a(Re) f + b(Re) f . (6) 

In this form the first term was said to represent the heat 
transfer through the laminar boundary layer on the front 
portion of the cylinder. The second component, according 
to Pand, represented the heat transfer in the separated 
region in the rear half. 
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Fand mentioned that in Ref. 18, Van der Hegge 
Zijnen proposed a similar form, and that Richardson in 
Ref. 19 stated that the term representing heat transfer in 

0.67 

the separated region is a "function of (Re)^ Fand sum- 
marized these developments and arrived at the following 
correlation formula, 

0.50 0.67 

(Nu) f = a(Re) f + b(Re) f . (7) 

He stated that Perkins and Leppert [20] correlated their 
data for water and ethylene glycol with formula (7). How- 
ever, Fand correlated his data successfully with the fol- 
lowing formula, similar in form to (7), 

(Nu) f = [0.35 + 0 . 3 4 ( Re + 0 . 15 (Re )’* ] (Pr (8) 

0.58 

where the term O.lSCRe)^ represents the heat transferred 
to the separated region. Fand opined that this formula was 
more "related to the physical processes involved" than 
McAdams' formula, and did, in fact, represent his data 
equally well. Fand further stipulated that equation (8) 
was only valid over the range of Reynolds numbers between 
200 and 100,000. 

c. W.H. Giedt 

Giedt [15] experimented with a constant heat 
flux cylinder. The relatively large cylinder (4 inch O.D.) 
that he used permitted the investigation of Reynolds num- 
bers from 70,800 to 219,000. The turbulence factor of the 
wind tunnel used by Giedt was given as 2.25 %• The 
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constant heat flux surface boundary was generated by elec- 
trically heating a thin Nichrome ribbon that was wrapped 
helically around the cylinder. The temperature rise of the 
cylinder surface above the free stream air temperature was 
measured with thermocouples that were welded to the inner 
surface of the Nichrome. The heating circuit was composed 
of five turns of ribbon, the center turn being the test 
section and the two on either side functioning as guard 
heaters. By maintaining constant power, and measuring the 
temperature increase, Giedt determined the heat transfer 
coefficient from equation (2) 

Q_ 

h = 4r- (9) 

The center of the hollow cylinder was filled 

/ 

with glass wool to reduce conduction losses, which were 
determined to be a maximum of two percent. A static-pres- 
sure tap was also installed on the cylinder. Then, by ro- 
tating the cylinder, both heat transfer and pressure 
information were collected at any desired angle to the flow 
direction. It was stated that in general, data were taken 
on one half of the cylinder circumference but that several 
measurements were made around the entire cylinder circum- 
ference with no substantial differences noted. 

Giedt performed a heat balance on a differential 
element of Nichrome ribbon and thereby accounted for heat 
loss by circumferential conduction through the Nichrome 
ribbon . 
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As Schmidt and Wenner [13] had found previously, 
Giedt again noted that heat transfer is a function of 
Reynolds number and "location along the cylinder circumfer- 
ence." The other major conclusions of Giedt are reiterated 
here for subsequent reference, as follows: 

1. Maximum values of h occur approximately ten degrees 
further back along the cylinder circumference than 
the minimum static-pressure values. 

2. Both minimum h values and pressures move back along 
the cylinder with increasing values of Re. 

3. For values of Re up to around 100,000, maximum h Q 
values occur at the 0 degree and 180 degree points. 

At values of Re = 1*J0,000 and up, h at the rear 

VJ 

stagnation point is greater than at the forward 
stagnation point. 

5. At values of Re from 1*10,000 to 219,000, distinct 
maximum h values occur between 110 and 115 degrees. 

6. For this same range of Re, two distinct minimum hg 
values occur. The lower of the two is between 80 and 
93 degrees on the front half of the cylinder. The 
second is between 1*10 and 150 degrees on the back 
half . 

7. Also for this same range of Re, inflection points 

occur in the pressure distribution curves at the same 

locations as the minimum h. values just noted. 

u 

It is noted that conclusion 3 above agrees 
identically with Schmidt and Wenner [ 1 3 1 • 
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Giedt made a significant comparison of his heat 

t 

transfer results with the skin-friction and static-pressure 
distributions given by Fage and Falkner [6,7]. In his dis- 
cussion for the runs with Re = 140,000 and 219,000 he 
stated that "the first minimum of the Nusselt number curves 
coincide in location with the first inflection points of 
the corresponding pressure curves. These inflection points, 
according to Fage and Falkner, mark the location of the 
minimum of the skin friction. The fact that the point unit 
heat transfer coefficients drop to a minimum and then 
abruptly rise near these inflection points also substanti- 
ates their conclusion that transition from laminar to tur- 
bulent flow in the boundary layer occurred at this same 
location." Giedt also concluded that "the second inflec- 
tion points which occur between 14,0 and 150 degrees in the 
pressure-distribution curves for these same values of Re, 
coincide with the locations of the second minimum hg val- 
ues." Giedt stated that based upon Fage and Falkners ' work 
"these second inflection points would also indicate the 
locations of second minimums or zero values of the skin 
friction." Giedt finally expressed a parallel between the 
location of the peak hg value in the vicinity of 110 degrees 
to the peak in the skin friction curve that rises after the 
first minimum. 

d. H.C. Perkins, Jr. and G. Leppert 

Perkins and Leppert experimented with two dif- 
ferent techniques using a constant heat flux cylinder. The 
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cylinders in both instances were small stainless steel 
tubes that were mounted in a test channel, heated uniformly 
by their electrical resistance to an applied, direct cur- 
rent, and cooled in a cross flow of fluid. 

Their first technique was especially intrigu- 
ing as it employed a constant heat flux cylinder to deter- 
mine the average heat transfer coefficient. As described 
in Ref. 20, a single thermocouple was located along the 
center-line of the cylinder, and the measured center-line 
temperature was mathematically shown to equal the mean 
temperature of the inside surface of the tube. This tern- 

P 

perature was used to calculate the temperature drop through 
the tube wall and thus to determine the surface mean temper- 
ature of the external surface of the tube. The average 
heat transfer coefficient was calculated from the latter 
surface mean temperature. The cooling fluid was water and 
ethylene glycol. 

Perkins and Leppert correlated the data from 
this first approach at Reynolds numbers from 40 to 100,000 
and Prandtl numbers from 1 to 300 using several methods of 
correlation. They were most satisfied with results ob- 
tained using the methods of Sieder and Tate for pipe flow 
and Vliet and Leppert for flow around spheres. The varia- 
tion in viscosity across the boundary layer was taken into 
account by means of a viscosity ratio factor. Other fluid 
properties were evaluated at the bulk temperature. 
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The other formulas used to correlate this data 
were of two forms, the one proposed by Richardson [19], 
equation (7), and the one proposed by Douglas and Churchill 
[17] , equation (6). Perkins and Leppert reported that both 
formulas, after being modified to include the Prandtl num- 
ber and the viscosity ratio, correlated their data well. 

It was stated that equation (6) as modified might be limited 
to Reynolds numbers not greater than 100,000. 

Perkins and Leppert also illustrated that the 
latter formulas after modification correlate the data of 
some others for a Reynolds number greater than 10. Fand 
[16], however, reported that equation (7) in its modified 
form predicts average Nusselt numbers that are 5 ^- 76 % 
higher than the experimental results that he obtained. 

Fand suggested that the "free stream turbulence, secondary- 
flows and boundary-layer effects associated with the rec- 
tangular duct in the apparatus of Perkins and Leppert are 
the cause of this discrepancy." 

In the second experimental technique, described 
in Ref. 21, Perkins and Leppert measured local heat trans- 
fer coefficients. The surface temperature measurement 
method and cooling fluid were different. This time a ther- 
mocouple was attached to the inside wall of the tube and 
utilized to measure temperature as a function of angle. 

The cooling fluid, water, had Prandtl numbers from 1 to 7 • 

A pressure pickup was also built into the cylinder to pro- 
vide a pressure measurement capability. Local heat 
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transfer coefficients were determined at Reynolds numbers 
from 2000 to 120,000. 

Both references reflect the significant effort 
made by the two experimenters with regard to blockage 
effects. Potential flow theory was utilized extensively 
for this analysis. 

An analysis was carried out [21] for predicting 
the heat transfer that occurs prior to separation using 
"integral methods with the Pohlhausen velocity profile in 
the momentum boundary layer and a third order temperature 
profile in the thermal boundary layer." 

Perkins and Leppert found [21] that "on the 
rear part of the tube, in the separated region, the tempera- 
ture underwent rapid and large fluctuations." An oscillo- 
graph was used to record "not only the temperature-time 
dependence at each angle in the separated region but also 
provided an accurate way to determine the time averaged 
temperature on that part of the tube." The local heat 
transfer coefficients were calculated using the temperature 
rise above the ambient air temperature, the measured, con- 
stant heat flux, and the application of Newton's law of 
cooling. The viscosity ratio correction was employed to 
develop local Nusselt number correlations. 

The turbulence level was reported to be about 
2.9 per cent. It was demonstrated, however, that the data 
collected using both techniques "appears more in line with 
a 1 per cent turbulence level." 
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e. Jerome R. Jarcy 

Jarcy conducted an experimental study of heat 
transfer from constant heat flux cylinders at relatively 
high temperatures. Reference 22 describes his study of flow 
field, drag and heat transfer characteristics of a mica 
cylinder. The cylinder was heated electrically using 
Nichrome ribbon wound around the cylinder. The average 
surface temperatures were from 250°P to 1100°F above the 
ambient air temperature. Jarcy reported that the flow for 
his experiment was subcritical and that separation was 
laminar. He determined this from the pressure distribution 
around the cylinder. Velocities involved were between 40 
and 160 feet per second. 

Jarcy’ s conclusions stated that "When a cylin- 
der is heated, the value of the dynamic and total pressures 
in the wake of the cylinder decreases as compared to a zero 
heating rate." 

f. P.D. Richardson 

Richardson considered [19] the separated flow 
region in the rear of a bluff body immersed in a transverse 
stream. He was particularly concerned with turbulent flows 
at high Reynolds numbers that do not reattach after separa- 
tion. Richardson concluded that heat transfer at the rear 
stagnation point "is proportional to the two-thirds power 
of Reynolds number." He also pointed out that it is known 
from the photographic record of Page and observations of 
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Roshko that in the separated flow region variation occurs 
in the direction of the axis of the cylinder. 

Richardson described an approximate method for 
calculating the convective heat transfer in the separated 
flow region of a bluff body and applied this method to a 
cylinder. This method accounts for the separate physical 
mechanisms of heat transfer that exist in the laminar 
boundary layer and separated regions. His method represents 
heat transfer in the separated region as being proportional 
to the two-thirds power of Reynolds number. He mentioned 
previous observations "that separated flow develops first 
at the rear stagnation point and then spreads forward round 
the cylinder at a decreasing rate with increase of Reynolds 
number until a position about 80 deg from the forward stag- 
nation point is reached." He reported that the initial 
development apparently occurs at a Reynolds number of 10. 

In Ref. 23 he stated that blockage effects are felt within 
the separated region. Richardson compared his results with 
the experimental data of others and reported good agreement, 
g. H.H. Sogin and V.S. Subramanian 

Reference 24 describes an experimental study of 
mass transfer from naphthalene cylinders in a cross flow of 
air by Sogin and Subramanian. They credited this technique 
to Froessling's work in 1938 and 1940, which made use of 
the analogy between constant-property incompressible heat 
transfer and low concentration incompressible mass transfer. 
The specific technique used in their investigation was said 
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to be the same as that of Christian who measured local 
rates of mass transfer from hollow napthalene cylinders in 
axisymmetric flow. Their measurements spanned a range of 
Reynolds numbers from 122,000 to 3^2,000. Their results 
before boundary layer separation were compared with the 
approximate analytical solutions of H. Schuh and H.J. Merk 
as well as other experimental works. 

2 . Some Available Theoretical Work 

Because of the complexity of the flow pattern once 
the laminar boundary layer separates or becomes turbulent , 
analytical solutions have not been found that predict local 
heat transfer over the entire cylinder. They generally are 
limited to approximate methods for calculating heat transfer 
in the front region of the cylinder extending from the for- 
ward stagnation line back to about eighty degrees. In 
Ref. 12, Spalding and Pun have compiled a summary, classi- 
fication, and comparison of fifteen theoretical analyses of 
heat transfer in this region for a constant temperature 
cylinder. Included are the analyses of Froessling, Seban, 
Squire, Eckert, Schuh and Lighthill. Spalding and Pun 
emphasized the need for experimental data collected at a 
low free stream turbulence level, to act as a standard of 
comparison and assist in the verification and evaluation of 
the proposed solutions. 

Seban and Chan present in Ref. 25 various methods 
for predicting heat transfer to laminar and turbulent 
boundary layers. The variable wall temperature case is 
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considered and, in particular, the constant heat flux case. 
Lighthill’s method for laminar boundary layers with almost 
constant properties is also included in the report and 
"recommended for general application, although in certain 
instances, the error in the predicted heat transfer may 
attain twenty per cent." 

Reference 26 contains Schuh' s method for calculat- 
ing laminar heat transfer on cylinders having either a con- 
stant wall temperature or a constant heat flux surface. 
Schuh compares the former with the data of Schmidt and 
Wenner [13] and the latter with the data of Giedt [15]. 

Good agreement was shown at Reynolds numbers below 100,000 
with deviations otherwise. Schuh reports that these devia- 
tions were pronounced at and near the forward stagnation 
point, with the experimental values up to 20# higher. 

In Ref. 27, Squire’s exact solution for forced con- 
vection heat transfer at the forward stagnation point of a 
cylinder is contained as an unpublished inclusion. His 
resulting value of Nusselt number is given as 



Nu = a. 



8 X D : 



v 



( 10 ) 



where a 3 is a function of the Prandtl number (a table of 
values of a 3 is provided), and is defined by 
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The velocity outside the boundary layer is given 



as 



u 



1 



4Ux 
D » 



( 12 ) 



where x is the distance from the forward stagnation point 
and U is the velocity of the uniform stream. Then, 






^ _ 4U 
x D 



(13) 



and taking for example, an average Prandtl number for air 
as 0.71, one obtains the result 



Nu 



0.995 



f UD ] % 

l v J 



(14) 



3. The Effects of Free Stream Turbulence on Heat 
Transfer Rate 

It has long been recognized that free stream turbu- 
lence has a significant effect on the heat transfer from a 
heated bluff body in cross flow. Some of the major efforts 
to resolve and account for these effects have been that of 
J. Kestin (Refs. 28, 29, 30, 31), P.P. Maeder [29], R.T. 
Wood- [31], R.A. Seban (Refs. 32, 33) and W.H. Giedt 
(Ref. 34). Kestin [28] describes a dual effect, one a 
broad change in the heat transfer rate resulting from the 
change in flow pattern, and the other a possible local ef- 
fect at a given angle. Some of his other conclusions are: 
a. The magnitude of the heat transfer variation due to 
the presence of free stream turbulence intensity 
seems to be larger at lower Reynolds numbers. 
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b. The turbulence intensity local effect seems to be 



an increase in the heat transfer rate. 

c. The primary local effect occurs in the heat transfer 
"across a laminar boundary layer developed under the 
influence of a favorable pressure gradient" and is 
"of the order of 10 per cent." For a laminar boundary 
layer near the forward stagnation line of a cylinder, 
the local effect is more pronounced and is "of the 
order of 70-80 per cent for an intensity of turbu- 
lence of 2. 5-3-0 per cent." 

d. The local effect on the heat transfer rate is non- 
existent across a turbulent boundary layer or for 
separated flow. 

Kestin mentions [30] that all laminar boundary 
layers "carry stochastic oscillations whose amplitude and 
frequency spectrum is related to the turbulence intensity 
in the free stream," and that if these disturbances do not 
amplify, the Tollmien-Schlichting theory of stability 
classifies the boundary layer as laminar. 

Giedt [3^] also found an increase of local heat 
transfer rates due to free stream turbulence, including a 
maximum 33-3 per cent increase over the calculated value at 
the forward stagnation point, when a rope net was placed 
about two feet upstream of the leading edge of the cylinder. 

Seban [32] summarized the effect of free stream 
turbulence as an over-all heat transfer increase due to 
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effects such as a lower Reynolds number for transition to 
turbulent flow which permits higher heat transfer rates, an 
increase of heat transfer rate to the laminar layer, and an 
alteration of the character of the flow in the separated 
region which may ultimately cause a movement of the separa- 
tion point. He also reported that for flow with a laminar 
boundary layer the local recovery factor doesn't vary due 
to free stream turbulence alterations. 

Regarding the nature of turbulence, Seban measured 
a longitudinal component intensity and spectral distribu- 
tion. Near the stagnation point he included a previous 
suggestion that a Goertler type instability might occur. 

More recently Kestin and Wood [31] found that the 
instability near the stagnation line of a blunt body is 
different from Goertler' s. As a result of this instability, 
Kestin and Wood reported that "the flow field acquires a 
system of evenly distributed, concentrated vortices arrayed 
in the axial direction and wrapped around the upstream por- 
tion of the body; their spacing is inversely proportional 

to the square root of the Reynolds number. The factor mul- 
-i- 

tiplying Re 2 depends on the turbulence intensity, Tu, of 
the free stream and has a nonvanishing value even for Tu=0, 
decreasing with an increase in Tu." Therefore these authors 
conclude that "the flow field in the laminar boundary layer 
outside a stagnation line is not two dimensional.... Thus 
the effect of free stream turbulence on heat and mass- 
transfer rates from the forward portion of blunt bodies is 
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indirect. First, a system of vortices is formed under all 
conditions, and this enhances the transfer rates. The 
latter are further enhanced when the intensity of turbu- 
lence is increased because the spacing is decreased 
thereby . " 

4 . Correlation Formulae For Forced Convection Heat 

Transfer 

Although there does not exist a complete analytical 
solution for heat transfer around a cylinder due to the 
complexity of the flow in the region where the boundary 
layer is turbulent or has separated, design requirements 
for this information must be satisfied. The answers are 
generally provided by means of empirically determined cor- 
relation formulae. All correlations are for average Nusselt 
numbers and not local Nusselt numbers. Although the present 
work is concerned with the latter, this discussion of cor- 
relations is included for the sake of completeness. The 
contributions in this branch of the study have been numer- 
ous and tailored for almost every conceivable combination of 
flow conditions. Perhaps the best known, and one of the 
earliest, was the work of Hilpert [35] in 1933 whose data 
included Reynolds numbers from about 500 to 200,000. Unfor- 
tunately, the turbulence intensity of his wind tunnel is 
unknown . 

Richardson in Ref. 36 concluded that Hilpert' s 
measurements "do not correspond well with present-day under- 
standing of the heat transfer distribution around a circular 
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cylinder." As mentioned before, Richardson subscribed to 
a correlation formula with a separate term to represent 
heat transfer in the separated region, that is proportional 
to Reynolds number to the two-thirds power. 

In Ref. 37 McAdams has put together a compilation 
of much of the work in this area. 

After a review of previous work. Van der Hegge 
Zijnen presented [18] a new correlation formula for combined 
free and forced convection heat transfer from cylinders to 
air . 



i 



ill 



Ill . THEORY 



The overall objective of this study was to experimental- 
ly determine local, convective heat transfer coefficients 
around a right circular cylinder placed in a cross flow of 
air exhibiting a low turbulence intensity (~0.5#)* To ac- 
complish this objective, an experiment was designed that 
embodied several basic theoretical concepts. These concepts 
include the governing equations for the flow of air around 
a two-dimensional cylinder mounted vertically in a closed, 
low-speed wind tunnel. Also involved are the relations 
that govern the heat generation technique employed and the 
modes of heat transfer to the environment. 

A. DESCRIPTION OF THE FLUID MOTION 

Figure 1 is a schematic diagram of the system consid- 
ered. Location one is upstream, location two is the 




Test Section 



Diffuser 



Figure 1. Schematic Diagram of the Air Flow System 
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entrance to the test section, and location three is the 
center of the test section. All pressure pickups that 
were mounted either in the wind tunnel or on the cylinder 
were within this system. The applicable equations for 
calculating the velocity in the test section (section 3) 
are (see nomenclature for definitions of terms): 



1. Continuity Equation (clear tunnel, without cylinder) 

PjAjUj = P 2 A 2 U 2 (15) 

2. Bernoulli’s Equation (clear tunnel, without cylinder) 




t i p u ! 

2 2 



( 16 )' 



3. Dynamic Pressure at Any Location, i 



q i - I <>i u i 



(17) 



4. Pressure Coefficient at Any Angular Location on the 
Cylinder, 0 



P - P 

c = IS li 

p %p 3 u i 



( 18 ) 



It was assumed that the density of air was constant 
throughout the system at any instant of time. This assump- 
tin implies that the flow was incompressible, and since 
density is a function of temperature, it also implies that 
any temperature rise through the system was uniform. 
Specifically, 



Pi = P 2 = P 3 = 
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P 



(19) 



Combining equations (15), ( 1 6 ) and (19) and noting that the 
Aerolab tunnel had a contraction ratio of A 2 /A 1 = 0.1, 
it is found that 



(Pj-P 2 ) = APi-r 0.99 



pu 



( 20 ) 



Ideal test section behaviour may be described as the 
state existing if the static pressure and velocity at loca- 
tion three are identical to those at location two. The 
static pressure variation across the test section was neg- 
ligible for the purposes of the experiment. Therefore the 
static pressure at location two, P 2 ,was assumed to be iden- 
tical to P 3 . 

The dynamic pressure at location three is given as 

,2 



q 3 = Japu, 



( 21 ) 



Ideally, U 3 = U 2 , and the ideal dynamic pressure at loca- 
tion three may be written as 



q ideal ^ pU 2 * 

Combining equations (20) and (22) 
tt 2 A Pi - 2 

^ pU 2 = 0.99 = q ideal 



( 22 ) 



(23) 



is written as a second result. The numerator can be meas- 
ured directly. 

It is now convenient to manipulate equation (18) into 
a more suitable form. Since the static pressure was 



assumed to be constant, the numerator may be rearranged as 
follows 

P 9 * P 3 * P 9 - P a • W 

This latter form can be measured directly and may be desig- 
nated as 



P - P = AP 
o 2 cylinder * 



(25) 



The denominator of equation (18) is sometimes described as 



^actual 



( 26 ) 



The ratio of q actual to q ideal is termed the tunnel cor- 
rection factor, X, and is a function of the particular wind 
tunnel. For the wind tunnel used, 

X = - a - ct - u - a - 1 = 0.99 . (27) 

q ideal 

Now, combining equations (18), (24), (25), (26) and (27) 
the pressure coefficient is found to equal 

C - . (28) 

P Ar x _ 2 

The free stream velocity is found by combining equations 
(23), (26) and (27). The final result is 

/ 2 ° A Pi _ 2 

U 3 =/ — . (29) 

3 P 

The Reynolds number may then be calculated directly, 

U 3 D 



4 5 



Re 



v 



(30) 



Since the above results were derived from equations 
one and two for a clear tunnel, blockage effects due to the 
presence of a bluff body in a constrained channel have not 
yet been included. Blockage corrections to the free stream 
velocity and Reynolds number are computed in Appendix A. 

B. THERMAL ENERGY GENERATION AND TRANSFER 

Thermal energy was generated by energizing thin Nichrome 
ribbon with a direct current, I. By means of the uniform 
electrical resistance of the Nichrome, R, electrical energy 
was transformed uniformly Into thermal energy. This tech- 
nique was described by Wood [1] as a Joulean heating rate, 
where 

Q = I 2 R. (3D 

Once generated, the thermal energy was dispersed by all 
three modes of heat transfer. Convection was the primary 
mode, and forced convection in particular was the mode of 
interest. References 38-40 support the assumption made 
that external free convection could be neglected. 

Convection internal to the cylinder was considered to 
be negligible, as the cylinder was filled with foam rubber 
which prevented any major circulation of air. 

It is helpful to consider an intuitive account of the 
external forced convection involved in the cooling of a 
heated cylinder immersed in a cross flow of air. At the 
forward stagnation point the cylinder receives the full 
strength of the wind and, as one would expect, is 



46 



relatively cooler. Proceeding back along the sides of the 
cylinder one would expect a continued decrease of the cool- 
ing effect of the wind and a gradual increase in the wall 
temperature of the cylinder. Once the boundary layer 
separates or becomes turbulent it seems logical that the 
cylinder wall would rapidly begin to be cooled again, by 
the turbulent packets of fluid "scrubbing" the surface and 
transporting thermal energy from the cylinder surface. 

The concept employed for quantitatively evaluating the 
external forced convection was the same as that used by 
Giedt [15]. It consists of applying Newton’s law of 
cooling. 



to the previously discussed expression, 

Q « I Z R . 

(See nomenclature for definitions of terms.) Then the basic 
expression for the local heat transfer coefficient can be 
written as 



Conduction of thermal energy out the ends of the cylin- 
der was prevented by the installation of guard heaters. 

Since data were collected at steady state conditions, radial 
conduction of heat energy into the interior of the cylinder 
was neglected. Circumferential conduction of thermal 




(32) 



h 



I 2 R 



(33) 
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energy through the- Nichrome was considered. This transfer 
should vary with position since the temperature gradient 
is a variable of position. Near the forward stagnation 
point, plus or minus about 20 deg, the temperature gradient 
should be small and circumferential conduction should have 
a minimal effect. However, one would expect that farther 
around the cylinder the trend would eventually be a flow of 
heat into a relatively cool region. A net flow of heat 
into an element in this region augments the heat generated 
within the element, causing the wall temperature to be 
higher than it would be otherwise. In the vicinity of the’ 
separation point the temperature gradient should be steep 
and the effect of circumferential conduction would be 
maximized. A net flow of heat out of an element in this 
region of relatively high temperature would tend to indicate 
a higher heat transfer rate than an ideal case without 
circumferential conduction. 

Radiation of electromagnetic thermal energy from the 
cylinder to the tunnel walls was taken into account and 
treated as a loss. Radiation losses were relatively small 
because of the small temperature differences involved, but 
were not considered to be negligible. 

An elemental heat balance has been developed for an 
element of Nichrome ribbon to derive the detailed expres- 
sion required for calculating the local heat transfer co- 
efficient, hg . The resulting' expression contains all 
applicable modes of heat transfer. The components of this 
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Figure 2. Heat Transfer for an Element of Nichrome Ribbon 

expression are depicted in Figure 2. (See nomenclature for 
definitions of terms.) The heat balance is first written 



Q k + Q i - Q h " Q rad " Q k " 33T Ax " °' 



( 34 ) 



It is noted that 



Q = A dT 

y k R n n 3x » 


(35) 
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(36) 


Q, =h fl A (T -T ) , 
h o s 0 0° * 


(37) 


Q , = ae (T^-T** ) A , 
rad 0 co y s * 


(38) 
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and 




Substituting into the heat balance one has 



( 39 ) 
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Rearranging 
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Rewriting in a finite difference formulation 



h 0" 



k A 
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(42) 



This formula is the core of the computer program in Appen- 
dix B that was written for data reduction. 

The above expression doesn't account for aerodynamic 
heating of the cylinder, and for the continuous rapid heat 
ing of the ambient air by the wind tunnel boundaries and 
impeller at velocities above 125 feet per second. Because 
the former effect tended to raise the unheated cylinder 
above its initial temperature and the latter caused the 
slower responding cylinder to lag the ambient air tempera- 
ture, the two effects labored in opposition. Their com- 
bined, resultant effect was calibrated out of the data as 
shown in Appendix C. 
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IV. EXPERIMENTAL APPARATUS 



A. LOW SPEED WIND TUNNEL 

The experimental study was carried out in the low speed 
Aerolab wind tunnel, located in the Aeronautics Laboratory 
in Halligan Hall, U.S. Navy Postgraduate School. A schemat 
ic diagram of the wind tunnel is provided in Figure 3. 

The turbulence intensity level of the wind tunnel was 
measured using a hot wire anemometer and its associated 
electronic equipment. The turbulence level ranged from 
0.67 for the low speed run to 0.^3 for the high speed runs. 
The details of these measurements and the equipment used 
are described in Appendix D. 

The wind tunnel is powered by a 100 hp electric motor, 
and has a four speed transmission. All experimental runs 
conducted at a Reynolds number less than 350,000 were in 
third gear and the high speed runs were conducted in fourth 
gear. The tunnel was capable of achieving a maximum speed 
of 200 mph with a clear test section. However, with the 
cylinder in place, the peak speed was only about 172 mph. 

Air speed information was determined by observing the 
pressure drop between two static pressure manifolds which 
indicated static pressure in the settling chamber (Figure 1 
location one) and near the entrance to the test section 
(Figure 1, location two). This static pressure difference, 
AP 1-2 , was read on the water-filled air speed manometer. 
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Corner vanes Drive Motor 






Figure 3. Schematic Diagram of the Wind Tunnel 
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Section B-B Section A-A 



The manifold for each location was connected to four static 
pressure pickups installed in the wind tunnel walls at that 
location. Setting the zero and reading the air speed 
manometer required a great deal of consistent, meticulous 
effort to obtain acceptable accuracy. 

B. CYLINDER DESIGN AND CONSTRUCTION 

All data for this project were collected with a 4 . 47 
in (o.d.) cylinder constructed from a cast acrylic resin 
tube. A second cylinder made from phenolic tubing has been 
designed and partially constructed, and is expected to be 
an improved model. 

As is apparent in the descriptions which follow, the 
design concept of both cylinders was the same as that of 
Giedt [ 15 ] with only some physical alterations. 

1 . The Cylinder Made From Acrylic Tubing 

The test cylinder was 32 in long and spanned the 
entire height of the wind tunnel test section. Its 4.47 
in (o.d.) diameter and 0.30 in wall thickness were the final 
dimensions after being machined into round from an original 
4.625 in (o.d.). 

The choice of a relatively large cylinder diameter 
permitted sufficiently large Reynolds numbers for examining 
all flow regimes in spite of a fairly limited air speed. 

A circuit for generating a constant heat flux on 
the surface of the cylinder was installed by helically 
wrapping Nichrome ribbon around the circumference of the 
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cylinder. The Nichrome ribbon was 0.375 in wide and 0.0031 
in thick, and had a resistance of 0.456 ohms per foot of 
length. This particular size was selected primarily on the 
basis of availability. It also possessed a sufficient total 
resistance to accommodate enough current to drive the sur- 
face temperature well above the ambient air temperature, 
even for high speed flows. 

To facilitate the wrapping process and to promote 
a smooth surface, the cylinder was grooved 2.5 threads per 
in, 0.375 in wide, and 0.005 in deep. This left a separa- 
tion of only about 0.025 in between individual wraps, giv- 
ing a close approximation to a uniform heat flux. The ends 
of the ribbon were led through slots from, and back into, the 
interior of the cylinder where they were connected to power 
leads. The ribbon was bonded to the cylinder with RTV-108, 
silicone rubber adhesive. 

The main circuit consisted of 39 wraps and had a 
total resistance of 20.527 ohms as measured on a Rosemount 
commutating bridge, model 920A. Variations in the resist- 
ance due to an increase in temperature were considered 
negligible over the temperature range involved. This cir- 
cuit was located approximately in the center of the tube, 
beginning 9*125 in from the base and extending up 15*125 
in of the tube's length. 

Two smaller guard heater circuits were similarly 
installed above and below the main circuit. These consisted 
of three wraps of Nichrome ribbon whose ends were led 
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through separate slots. Each pair of guard leads were 
then connected to a separate pair of power leads. This 
arrangement allowed the guard heater circuits to be con- 
trolled individually. As shown in Figure 4, the guard 
heater circuits were laid down in the same continuous 
groove as the main circuit and adjacent circuit slots were 
separated by only 0.5 in. The smaller circuits each had a 
resistance of about two ohms. 

While the guard heater circuits were the primary 
means of precluding heat conduction losses out the ends of 
the tube, the wall thickness at each end was also shaved to 
0.0625 in to minimize conduction losses. 

Temperature information was measured using premium 
grade copper-constantan bead thermocouples made from 0.005 
in wire . One thermocouple was shielded and located on a 
screen in the settling chamber, where the velocity was a 
tenth of the test section velocity. It therefore provided 
a true measure of the ambient air temperature. It was 
referenced against an ice junction, and acted as a refer- 
ence for the remaining eight thermocouples which were in- 
stalled on the cylinder surface in the following manner. 

They were led through drilled holes in the cylinder wall 
and were spot welded to the inner surface of the ribbon just 
prior to the sealing of the ribbon to the surface. These 
thermocouples directly indicated the temperature difference 
between the cylinder surface and the ambient air. Their 
location on the cylinder is pictured in Figure 4. 
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Figure 4. Test Cylinder 
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The test section was designated as the center wrap 
of the main circuit. It consisted of four thermocouples 
attached at 90 deg intervals. The remaining four thermo- 
couples were used solely for guard heater control. Two 
were located on the center wrap of the two guard circuits. 
The other two were located near the edges of the main cir- 
cuit and provided the comparative temperature information 
required for the guard heater system to function. 

To avoid any contamination or damage to the test 
cylinder, the actual test thermocouples were not calibrated. 
Instead, an identical thermocouple was welded to a small' 
strip of Nichrome. This representative test thermocouple 
was then calibrated in the Rosemount constant temperature 
bath system which had a platinum resistance thermometer 
connected across the Rosemount commutating bridge as a 
standard. An ice bath was used for the reference junction. 
Since the results showed close agreement with the published 
table for copper-constantan thermocouples over the applica- 
ble temperature range, the table was used for data 
reduction . 

An additional check was performed on a similar 
representative test thermocouple to insure that no extrane- 
ous EMF, as described by Davenport and others in Ref. 41, 
was generated due to its attachment to a direct-current 
heater. The test consisted of pulsing the small strip of 
Nichrome with direct current and observing the response of 
the attached thermocouple on a Honeywell chart recorder. 
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Only an exponential rise and fall of thermally induced EMP 
were observed. No step changes, which would be indicative 
of an extraneous electrical EMF, were observed. This con- 
clusion was in agreement with Davenport [ 4 1 ] for a bead 
thermocouple making point contact with an electrically 
heated substrate. 

Another simple experiment was performed to gain 
insight into the relatively larger thermal expansion of the 
acrylic tube in comparison to the Nichrome ribbon. The 
objective was also to determine an upper limit for surface 
temperature, based upon the endurance of the acrylic to 
resist a combination of continuous high temperatures and 
potentially deforming hoop stress which was applied by the 
cinching action of the ribbon. A third objective of the 
experiment was to determine the ability of the R.T.V. 
sealant to resist high temperatures. The details of the 
experiment are given in Appendix E. The results are brief- 
ly summarized as follows. Deformation is probable in still 
air at temperatures above 130°P, and the R.T.V. held solidly 
at- temperatures up to l80°F. 

As the R.T.V. cured after the Nichrome ribbon was 
bonded to the final test cylinder, a contraction of the 
sealant occurred, resulting in various peaks and recesses 
in the ribbon and therefore a somewhat lumpy surface rather 
than the desired smooth one. The effects of these flaws 
are discussed later in the thesis. 
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Two static pressure pickups were installed in the 



cylinder wall 0.75 in. below the lower guard heater and 
at 90 deg intervals. These were made of 0.040 in (o.d.) 
stainless steel tubing, and had an orifice of 0.032 in 
(i.d.). In this way, surface pressure was measured and 
input to one side of a U-tube water-filled manometer. The 
manometer was referenced against the static pressure at 
the entrance to the test section, P . 

After the thermocouples and pressure taps were in 
place, the cylinder was filled with Silastic, S-5370, RTV 
foam, a low density resilient silicone rubber. The foam • 
rubber was easily poured into the cylinder and quickly 
expanded, uniformly filling every crevice. It provided 
support for the fine thermocouple wires. Internal free 
convection was practically eliminated, and the foam was 
an excellent insulator, exhibiting a thermal conductivity 
of 0.026 BTU/hr-ft-°F. Unfortunately, a serious disadvan- 
tage of the foam was its large thermal capacity which 
considerably increased the delay before reaching steady 
state conditions. Long experimental runs were the out- 
come of this objectionable feature. 

The cylinder was mounted on a turntable base 
that could be rotated from outside the wind tunnel. A 
disk, that was marked every five degrees, rotated with the 
cylinder base. Therefore, temperature difference and 
surface pressure could be determined at any desired angular 
location . 
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Figures 5 and 6 show the cylinder mounted in the 



test section. 

2 . A Proposed Phenolic Test Cylinder 

A second cylinder made from 4 in (o.d.) phenolic 
tubing has been designed and partially constructed. It is 
based upon the same experimental concept as the first 
cylinder, but has several different and potentially better 
features . 

The coefficient of thermal expansion is not as 
large as the acrylic tube, but is closer to that of the 
Nichrome . Therefore, it is possible that a sealant would 
not be required over the entire length of the ribbon, if 
required at all. This would preclude the development of 
surface irregularities. 

The phenolic tubing can withstand much higher 
temperatures for long periods of time without experiencing 
deformation . 

A space for four additional thermocouples has been 
provided on one side of the cylinder in the test section 
region, resulting in thermocouples spaced every 30 deg. 

This feature will expedite data collection. 

At the time of this writing, the phenolic tubing 
has been machined and grooved to accept the three circuits. 
The rotating base has been modified to accept the smaller 
diameter tube, and electrical design requirements have been 
completed . 
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Figure 5. Test Cylinder Mounted In the Wind Tunnel Test Section 






Figure 6. Closeup View of Test Cylinder 
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C. ELECTRICAL CIRCUITRY AND EQUIPMENT 

All power leads and thermocouple wires were led down 
the interior of the cylinder and out of the wind tunnel 
through a hole in the cylinder base. Each pair of power 
leads for the three circuits was connected to individual 
Lambda model LK 3^5A FM 60-volt regulated power supplies. 

A Non-Linear Systems Inc., Model 48l digital voltmeter was 
used to monitor the power settings for any of the three 
circuits to the nearest millivolt. 

The thermocouple wires were connected to a ten posi- 
tion rotary switch: eight positions for the thermocouples 
on the cylinder, one position for the thermocouple in the 
settling chamber and the tenth position was used to include 
the ice junction. The output leads from the switch were 
connected to either a Honeywell chart recorder during the 
transient, warming up period, or to a Leeds and Northrup 
millivolt potentiometer for the accurate data 'collection 
at steady state conditions. 

Figures 7 and 8 illustrate the arrangement of the 
equipment . 
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Figure 7 . Arrangement of the Experimental Apparatus 
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Figure 8. Close Up View of the Electrical Equipment 



V. EXPERIMENTAL PROCEDURE 



In the development of a precise experimental technique 
for the test cylinder constructed for this study and others 
that might follow, there were three distinct phases of ex- 
perimentation, each spanning several months of research. 

Initially, the heat transfer characteristics of the 
cylinder were investigated at a single Reynolds number of 
153,000. During this time a large amount of preliminary 
data were recorded to determine if the temperature distrib- 
ution over one side of the cylinder was the same as on the 
opposite side so that data could be collected from one side 
only; to ascertain which side of the cylinder was prefer- 
able and would provide the smoother data, particularly in 
view of the surface irregularities mentioned previously; 
and to resolve a preferable pattern of rotation in conjunc- 
tion with the preferred side. 

The second phase of experimenting consisted of an in- 
vestigation of eight different Reynolds numbers, and was 
originally intended to provide the final results. This 
would, in fact, have been the case had the author been 
more experienced in the art of wind-tunnel testing. How- 
ever, at various times in the process of collecting this 
data, flaws in the experimental technique were discovered 
which detracted from the credibility of the results. For 
this reason, the author elected to perform an additional 
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series of runs which incorporated a detailed experimental 
procedure. The following tasks, which compose this proce- 
dure, are discussed in their approximate order of 
accomplishment . 

The air speed and U-tube manometers were carefully set 
to zero, and all pressure line connections were checked 
for tightness. 

The Nichrome ribbon wrappings on the test cylinder 
were cleaned before each run with acetone to remove any 
film of dust, oil or oxidation. The cylinder was oriented 
into its initial position and the accesses to the test 
section were secured. 

All electronic equipment was turned on to begin warming 
up. The initial voltmeter reading was recorded for the 
main circuit with no current flowing. This reading was 
invariably -O.O65 volts. 

The wind tunnel was put in the proper gear, turned on 
and the air speed increased to only about ten mph. This 
circulated and brought the air to a uniform temperature. 

The zero was set on the millivolt potentiometer and 
the initial readings for all thermocouples were recorded 
including the initial air temperature. Just prior to 
recording these and all subsequent temperature data, a 
check was made to insure that a proper water-ice interface 
existed at the ice junction. 

The air speed was increased to provide the Reynolds 
number desired for the run. Choosing the correct air speed 
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for this purpose was not obvious. The air temperature in 
the tunnel generally increased with time to some level, 
at which the time rate of temperature rise was a minimum. 

For higher velocities the air temperature increased to 
higher levels. Since the viscosity and density of the air 
are dependent upon the air temperature, prior knowledge of 
the resultant level of air temperature was required for 
setting the air speed of the run. 

The output leads from the thermocouple rotary switch 
were removed from the potentiometer and connected to the 
chart recorder. The voltage to the main circuit was then 
increased gradually to the desired power level for the run. 
This level was chosen to be the amount of power that would 
maximize the temperature increase of the cylinder surface 
above the ambient air temperature without surpassing the 
high temperature endurance limit of the acrylic tube. 
Maximizing delta T enhanced the degree of accuracy of the 
heat transfer data and tended to alleviate many of the 
surface irregularities by driving more thermal energy into 
the sealant, thereby expanding some of the recesses. 

Power was increased to each guard heater circuit until 
an approximately equivalent amount of current was flowing 
as in the main circuit; the voltage was about an order of 
magnitude less than the main circuit voltage. As the main 
circuit settled to the desired level of power, fine adjust- 
ments were made to the guard heaters. Their surface temper- 
ature should be maintained equal to the surface temperature 
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at the edges of the main circuit by comparing the guard 
heater thermocouple readings to those on the edges of the 
main circuit. 

After the thermal capacity of the acrylic tube and its 
foam rubber filler were saturated, a steady state condition 
existed. Then, all of the heat added was removed by convec- 
tion, less any losses due to thermal radiation to the tun- 
nel walls or circumferential conduction through the 
Nichrome ribbon. The progress of the surface temperatures 
toward their steady state values was observed and tracked 
on the chart recorder. It generally required about two 
hours to reach steady state, depending on the velocity of 
the air. The higher speed runs required up to three hours 
and the lower speed runs were sometimes steady after 1.5 
hours. 

The experiment was modeled on TRUMP, a computer program 
developed at Lawrence Radiation Laboratory, Livermore, 
Calif., for solving transient and steady-state temperature 
distributions in multidimensional systems. Reference 43 
is the manual for applying TRUMP. The time required to 
reach steady state was predicted reasonably well using the 
TRUMP model. 

After the cylinder was in a steady state condition, the 
thermocouple rotary switch output leads were again connect- 
ed to the potentiometer. The potentiometer was again 
mechanically and electrically set to zero, as it was just 
prior to any recording of temperature data, and the steady 
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state values of the surface temperature Increases above 
the ambient air were recorded for the initial position of 
the cylinder. The ambient air temperature was also record- 
ed. The cylinder surface pressure was measured by the U- 
tube manometers and recorded. 

The cylinder was then rotated five degrees to a new 
position. The time required to reach steady state at the 
new position was generally 15 to 30 minutes. This was more 
than adequate time to check the ice junction, collect the 
pressure data at this new position and zero the potentio- 
meter. Once steady state was reached at the new position, 
the temperature data were recorded as before, and the 
cylinder again rotated five degrees. 

In this way, temperature and pressure data could be 
obtained at any desired angular location. All data were 
collected on one side of the cylinder only. The experimen- 
tal runs required 8-10 hours each due to the large thermal 
capacity of the cylinder, and the small number of test 
thermocouples installed. 
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VI. PRESENTATION OF RESULTS 



Raw data, including temperature information in milli- 
volts and pressure information in cm or in of water, were 
read directly into the computer program in Appendix B. 

The density, thermal conductivity and kinematic viscosity 
of air were evaluated at the film temperature, T^, which 
was computed by 

(T _) = (T , ) + 0 . 5 ( AT) (43) 

f run amb avg avg 

for each run, where (T , ) is the average ambient tem- 
* amb avg & 

perature and (AT) is the average local surface tempera- 
ture increase at steady state. 



A. TABULAR RESULTS 

The computer program provided detailed results for all 
thirteen experimental runs. These are presented in tabular 
form on pages 72-84. The pressure coefficient, heat trans- 
fer coefficient, Nusselt number and Froessling number are 
provided at five deg intervals from the forward to the rear 
stagnation points for all runs except three runs which 
omitted the pressure information. Two of these runs are 
duplicates of others which include pressure information. 
Therefore, a complete set of results is presented at ten 
different Reynolds numbers. 

"DELTAX" is the arc length, "ASURF" is the surface area, 
and "QGEN" is the heat generation rate of a 5 deg element 
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FORCEO CONVECTION HEAT TRANSFER FROM A UNIFORMLY 
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FORCED CONVECTION HEAT TRANSFER FROM A UNIFORMLY 
HEATED CYLINDER IN A CROSSFLOW OF AIR 
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of Nichrome ribbon. "TU" is the turbulence intensity level 
at the given free stream velocity. "DELTA-T CORRECTION" 
is the correction for the resultant effect of viscous heat- 
ing of the cylinder and heating of the ambient air during 
high speed runs (Appendix C). 

B. DISCUSSION 

To illustrate the discussion, the Nusselt number and 
the pressure coefficient have been plotted as a function of 
their angular location in Figures 9 through 14. 

The results in general only confirmed the trend of 
those published by Giedt C 15 J with regard to approximate 
locations of minimum and maximum points and the angular 
pattern of heat transfer characteristics around the cylin- 
der. The magnitudes of the heat transfer rates are some- 
what less than those of Giedt. This is most probably 
explained by the lower turbulence level of the wind tunnel 
employed in this experiment ( ~0.5 % )• Also, the alignment 
of the cylinder was about one deg in error, as it leaned 
toward the rear of the test section. This discrepancy 
probably decreased the normal component of the wind strik- 
ing the cylinder, and may have consistently decreased the 
heat transfer rates. 

At the lower Reynolds numbers of 57,000 and 82,000, the 
heat transfer results were in close agreement with the 
resulting pressure coefficient distribution around the 
cylinder. The boundary layer remained laminar and 
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Nusselt Number 



o Re = 57 ,000 
□ Re = 223,000 
ARe = 350,000 




Figure 9* Heat Transfer Results at Reynolds Numbers of 
57,000, 223,000 and 350,000 
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Nusselt Numbers 




Figure 10. Heat Transfer Results at Reynolds Numbers of 
82,000, 261,000 and 406,000 
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Figure 11. Heat Transfer Results at Reynolds Numbers of 
153 , 000 , 309,000 and 495,000 
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Figure 12. Surface Pressure Distribution at Reynolds Numbers of 57,000, 223,000 
and 309,000 
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Figure 13. Surface Pressure Distribution at Reynolds Numbers of 82,000, 261,000 and 
406,000 



153,000 

495,000 
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Figure 14. Surface Pressure Distribution at Reynolds Numbers of 153,000 and 495,000 



separated at about 82.5 and 80 deg, respectively. This 
laminar flow regime was reflected in the heat transfer 
results by a single minimum value of the Nusselt number, 
which was located at the separation point. Also, maximum 
heat transfer occurred at the stagnation points. The cor- 
responding pressure distributions also depicted laminar 
flow. These were characterized by a minimum at 65 deg 
followed by a fairly small recovery in the form of a 
smooth, continuous rise, and a flattening out at the sepa- 
ration points. Additionally, there were distinct peaks in 
the two pressure coefficient curves at 85 and 80 deg, res-, 
pectively, just as the curve flattened out. 

Unfortunately, at the intermediate range of Reynolds 
numbers studied (82 ,000<Re<300,000) , the surface pressure 
distribution did not reflect the same conditions of flow 
inferred by the heat transfer results. As before, the 
pressure coefficient distributions were representative of 
subcritical flow, and indicated a laminar boundary layer 
separation at about 80-85 deg. The heat transfer results, 
however, depicted a region of transitional flow in the 
boundary layer, and therefore were representative of crit- 
ical flow. There were, for example, two minimum values of 
Nusselt number. The first minimum occurred at an angular 
location of about 85 deg, and marked the beginning of 
transition to a turbulent boundary layer. At this point 
the boundary layer may have separated, become turbulent, 
and then reattached at a subsequent point, having become 
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fully turbulent. This transition region, known as the 
"separation bubble", is exhibited by the angular distribu- 
tion of the local Nusselt number, at this range of Reynolds 
numbers. The second heat transfer minimum occurred at about 
130 deg, and represented the separation of the turbulent 
boundary layer. At this range of Reynolds numbers, and 
between the two minima described, maximum heat transfer 
occurred at an angular position corresponding to the re- 
attachment of a fully turbulent boundary layer. 

The lack of agreement of the pressure information can 
most probably be attributed to the location of the pressure 
pickups on the smooth acrylic tube and not in the region of 
the test section. 

At the higher Reynolds numbers of * 106,000 and * 195 ,000 
there was close agreement between the heat transfer and 
surface pressure results. Both indicated critical flow 
conditions. For each run the heat transfer results exhib- 
ited two minimum values of the Nusselt number. The first 
minimum values where transitional flow began were located 
at 90 and 82.5 deg, respectively. The second minimum 
values where the boundary layer separated were located at 
130 and 135 deg, respectively for the two runs. The points 
of maximum heat transfer where the boundary layers reat- 
tached and became fully turbulent were located at 100 and 
95 deg, respectively. 

For these same two Reynolds numbers the surface pres- 
sure distributions had minimum values at 80 deg, which 
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differed from all other experimental runs which had pressure 
minimums at 65 deg. Also, the rise of the pressure curve, 
from the minimum, was not smooth, but was characterized 
by a double inflection, which depicted the "separation 
bubble." The first inflection point represented the initial 
detachment of the laminar boundary layer. Its angular loca- 
tions of 90 deg and 85 deg for the two runs corresponded 
identically with the locations of the first minimum values 
of the Nusselt number. The second inflection point indi- 
cates reattachment and complete transition to a turbulent 
boundary layer. Its angular location of 100 deg in both 
runs, corresponded closely with the location of the maximum 
values of the Nusselt number. Also the pressure distribu- 
tions for these Reynolds numbers flattened out at angular- 
locations of about 130 and 120 deg, respectively, corres- 
ponding reasonably well with the second minimum values of 
the Nusselt number. 

In contrast to Giedt’s findings [15] for Reynolds 
numbers above 153*000, the value of hg at the rear stagna- 
tion "point was not always greater than its value at the 
forward stagnation point. 

It was not found that the angular location of minimum 
values of hg invariably increased as the Reynolds number 
increased. Instead, there were variations in the flow con- 
ditions from one run to the next. A higher velocity run 
at times produced an apparently smoother flow than one of 
the lower velocity runs. For example, the initial minimum 
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values of h q for Reynolds numbers of 153,000; 173,000; and 

223.000 were 90,100 and 100 deg, respectively. However, 
at a Reynolds number of 260,000, the first heat transfer 
minimum occurred at 82.5 deg. In the same way, the initial 
minimum value of hg for Reynolds numbers of 309,000; 

350.000 and 406,000 were 85 , 95 and 90 deg respectively. 

But apparently a very smooth flow existed at a Reynolds 
number of 495,000, since the first heat transfer minimum 
then occurred at 85 deg. This inconsistancy was probably 
caused by the surface defects of the test cylinder. These 
defects were less pronounced during the two high speed 
runs because of the higher surface temperatures involved. 

It is thought that the large thermal capacity of the 
foam rubber and acrylic in the interior of the tube fos- 
tered complicated internal flows of heat. These heat flows 
depended upon the magnitude and direction of the angular 
rotation of the cylinder from a previous steady state posi- 
tion to the present position. Also, the magnitudes of the 
steady state surface temperature differences varied, depend- 
ing upon whether the cylinder was rotated to a position 
from a previous steady state position, or was initially 
warmed up there and heated to a steady state condition. 
Specifically, and for example, rotating a test thermo- 
couple at five deg intervals gradually around a distance 
of 20 deg to a critical location seemed to amplify these 
temperature differences than if the cylinder had initially 
been warmed up with a test thermocouple in that critical 
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location. This was especially true when approaching and 
arriving at the first minimum value of the Nusselt number. 
Variations near the front of the cylinder were only about 
0.5°F or less, but were on the order of five °P at the 
critical position, and about one °F at the rear stagnation 
point . 

An uncertainty analysis was performed and is described 
in Appendix P. In this analysis the uncertainty of the 
Nusselt number was found to be about three per cent at the 
forward stagnation point, about six per cent at the rear 
stagnation point, and on the order of 25# at the critical 
point where the boundary layer either begins its transition 
to turbulent flow or separates. 

C. COMPARISON OP EXPERIMENTAL RESULTS WITH EXISTING THEORY 

The heat transfer at the forward stagnation point can 
be conveniently compared to Squire's solution, which gener- 
ates a Proessling number (P = Nu//Re) of 0.995 for this 
application. The Proessling number calculated from the 
present set of experimental results was about 0.93-.92 at 
Reynolds numbers from 153,000 to 261,000. At Reynolds 
numbers below and above this range, the Proessling number 
decreased to as low as .84, excluding the run at a Reynolds 
number of 495,000. This latter run had a surprising, low 
Proessling number of 0.76 at the forward stagnation point, 
which could not be explained. For this run the power 
source for the main circuit and the wind tunnel were both 
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in their maximum positions, which possibly placed an unac- 
ceptable strain oh the experiment. It is also curious to 
note that the Nusselt number at the forward and rear stag- 
nation points are identical for the two experimental runs 
of 406,000 and 495,000. 

For an evaluation of the data beyond the forward stag- 
nation point an appropriate theoretical solution was re- 
quired. All existing theoretical solutions only predict 
the heat transfer from the forward half of a cylinder, 
prior to the point of laminar boundary layer separation at 
about 80 deg. Spalding and Pun [12] have summarized fif- 
teen approximate solution methods. However, since all fif- 
teen are applicable to isothermal cylinders, none would 
provide a meaningful comparison for the data collected in 
this experiment with a constant heat flux cylinder. 

Spalding and Pun compared the various methods with that of 
Froessling who used an "exact", series method of solution 
limited by the number of terms employed. It was shown that 
Schuh* s method of solution for an isothermal cylinder 
ranked in the highest class of accuracy and within 1-3 per 
cent of Froessling* s solution. Schuh also published [26] 
one of the few approximate methods of solution for a con- 
stant heat flux cylinder. 

' In Fig. 15 the two theoretical solutions of Schuh are 
plotted to illustrate the significant difference in the 
heat transfer profile between the two types of cylinders 
usually employed. As one would expect, the two solutions 
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Nu//Re 




Figure 15. Theoretical Solutions of H. Schuh for Both 
Isothermal and Constant Heat Flux Cylinders 
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converge at the forward stagnation point, for in that vicin- 
ity the constant heat flux cylinder also has a constant wall 
temperature. Moving farther around the cylinder the ratio 
of relatively greater heat transfer from the constant heat 
flux cylinder continually increases. This is physically 
represented by an increasing surface temperature, which is 
required to provide an increasing temperature difference 
across the thermal boundary layer. The increased tempera- 
ture difference is required to drive a constant flux of 
thermal energy through the thickening thermal boundary 
layer. 

Some of the experimental results of Giedt [15] for 
heat transfer from a constant heat flux cylinder are com- 
pared to Schuh's solution in Pig. 16 . The nature of the 
divergence of Giedt' s results from the theoretical curve 
is most probably explained by the significantly higher level 
of free stream turbulence ( 2 . 25 % compared to 0.5%). 

In Fig. 17 some of the present results are compared 
with Schuh' s 'solution, and most of them agree within about 
six per cent. This agreement is largely due to the low 
turbulence level of the wind tunnel used (0.5%). These 
results indicate that slightly less heat transfer occurs 
than that predicted by Schuh. The results are closely 
grouped, implying that any experimental errors are fairly 
consistent . 
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Figure 16. Comparison of Giedt's Experimental Results to 
the Theoretical Solution of Schuh 
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Nu//Re 



O Re = 82,000 
□ Re = 153,000 
A Re = 223,000 
X Re = 259,000 
-f Re = 308,500 




Figure 17* Comparison of Experimental Results of the 
Present Investigation to the Theoretical 
Solution of Schuh 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



Heat transfer and pressure coefficients were determined 
around a heated, right circular cylinder placed in a cross 
flow of air. The heat transfer results were generally in 
agreement with the pattern of heat transfer found previ- 
ously for subcritical and critical flow. 

The heat transfer results were compared with the clas- 
sic results of Giedt [15], and indicated relatively lower 
heat transfer rates. This is most probably explained by • 
the considerably lower free stream turbulence level of the 
wind tunnel employed in this experiment. It has been pre- 
viously shown by other researchers that higher free stream 
turbulence levels markedly increase the heat transfer 
rates . 

The heat transfer results were also compared and showed 
good agreement with the approximate analytical solution of 

I 

Schuh for the laminar flow region on the front portion of 
the. cylinder . 

The surface pressure distribution did not in all ex- 
perimental runs support the heat transfer results. Specif- 
ically, in several runs the pressure information failed to 
indicate the "separation bubble" which represents the trans- 
ition in the boundary layer to turbulent flow, and depicted 
subcritical flow instead. In these instances the heat 
transfer results clearly represented critical flow effects. 
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It is possible that the pressure measuring system was not 
adequately sensitive to detect all of the surface pressure 
variations. However, it is more probable that the loca- 
tion of the surface pressure pickup was too distant from 
the test section. It is recommended that on the next cyl- 
inder constructed, the pressure taps be located in the test 
region . 

An uncertainty analysis was performed and showed that 
the uncertainty in the heat transfer results was about 
three per cent at the forward stagnation point, six per 
cent at the rear stagnation point and on the order of 2 5$ 
at the peaks in the heat transfer curves, where fluctuations 
seemed to be amplified. 

The heat transfer results were possible adversely 
affected by some defects in the test .cylinder, which can 
be avoided in a second cylinder being constructed. 

The easily deformable acrylic tube and its large ther- 
mal expansion coefficient as compared to Nichrome ribbon, 
and the curing properties of the sealant used to bond the 
Nichrome to the surface resulted in a number of surface 
irregularities. A plan for avoiding this defect is being 
implemented in the construction of the second cylinder. In 
the "Discussion" by B.H. Sage appended to Ref. 32, small 
surface irregularities are reported to have caused unusual- 
ly large variations in the heat transfer, especially in 
the rear portion of a cylinder. Seban [32] stated that 
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similar surface defects "disturbed the flow near the separa- 
tion line." 

It is recommended that the TRUMP computer program [43] 
be utilized to evaluate various insulating materials as a 
suitable replacement for the foam rubber. It should be 
possible to eliminate the complicated thermal flows and 
rotational behaviour, long transient warm-up period, and 
possible internal free convection as the acrylic tube 
usually expanded clear of the foam rubber filler. 

It is recommended that the cylinder vertical alignment 
be carefully checked before all experimental runs. The 
base plate moving parts and metal interfaces should be 
overhauled, sanded smooth and lubricated. The rotating 
disc should be marked to indicate every one deg; it is 
the author's opinion that in a tunnel with a low turbulence 
level, a sufficiently accurate study is possible that 
warrants this refinement. The wooden floor has become 
warped from heat and should be renewed. 

A major source of concern with regard to accurately 
measuring the free stream velocity was the air speed mano- 
meter. It is recommended that a hot wire anemometer be 
installed near the cylinder in such a manner as not to 
disturb the flow. This also precludes the problem of in- 
volving corrections for wind tunnel boundary constraints. 

It is also recommended that a separate study of the 
same cylinder be made using cholesteric liquid crystals as 
a temperature sensor on the surface of the cylinder. 
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These crystals indicate temperature through a change in 
color, and offer the possibility of providing a highly 
desirable visual display of the temperature field much as 
smoke streaks provide a visual display of various flow 
phenomena. 
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APPENDIX A 



Wind Tunnel Blockage Corrections 



The test section velocity obtained from the air speed 
manometer and its associated pressure pickups refers to a 
clear test section. In Ref. 42, Pope discusses the effects 
of the tunnel boundary restrictions in combination with a 
two-dimensional model mounted in the test section. The 
flow area is reduced by the physical presence of the model 
and its relatively slower moving wake. Both of these dis- 
turbances cause the air to accelerate as it flows past the 
model in order to preserve continuity. Hence the measured 
velocity must be corrected for solid blocking and wake 
blocking, respectively. 

Calculation of the solid blockage correction : 

This boundary correction was calculated using the 
equation of Pope [42]. 



e SB 

e SB 




n 2 




e SB = 0.008215 



(44) 



Calculation of the wake blockage correction : 

This boundary correction is a function of the drag 
coefficient, and had to be evaluated for each run. The 
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equation used is from Pope [42], 
were taken from Fig. 1.4 of Ref. 






e WB To 




Run 


Uncorrected 
Reynolds Nr. 


C D 


1 


54,785 


1.20 


2 


79,296 


1.20 


3 


97,522 


1.20 


4 


147,675 


1.20 


5 


166,381 


1.20 


6 


215,044 


1.20 


7 


215,185 


1.20 


8 


250,655 


1.05 


9 


251,928 


1.05 


10 


298,958 


0.95 


11 


340,095- 


0.85 


12 


395,887 


0.70 


13 


487,288 


0.30 



and the drag coefficients were 

5. 

(45) 



e WB 


e t = e SB 


0.03 


0.038 


0.03 


0.038 


0.03 


0.038 


0.03 


0.038 


0.03 


0.038 


0.03 


0.038 


0.03 


0.038 


0.026 


0.034 


0.026 


0.034 


0.023 


0.031 


0.021 


0.029 


0.017 


0.025 


0.007 


0.015 



Application of the total blockage correction 

The right column of the table provides the total block- 
age correction required for this experiment. Only the free 
stream velocity and the Reynolds number were corrected for 
boundary effects. The boundary corrections did not enter 
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into calculation of the pressure coefficients. 
Specifically , 

U corrected = U uncorrected (l+e^) 

Re corrected = Re uncorrected (1+e^) 
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APPENDIX B 



Computer Program for Data Reduction 



THIS PROGRAM INCORPORATES A FINITE DIFFERENCE SCHEME 
TO COMPUTE THE THERMAL ENERGY BALANCE AT EACH ELEMENT 
ASSOCIATED WITH A NODE OR DATA POINT. A NODE/DATA 
POINT LIES IN THE CENTER OF AN ELEMENT. AN INTERFACE 
BETWEEN ELEMENTS IS HALFWAY BETWEEN ANY TWO NODES. 
EACH ELEMENT IS 0.375 INCHES WIDE AND 0.003L INCHES 
THICK. THE ARC LENGTH OF AN ELEMENT IS EQUAL TO THE 
DISTANCE BETWEEN INTERFACES WITH ADJACENT ELEMENTS. 



UNITS EMPLOYED IN THE PROGRAM ARE: 

LENGTH-FEET 
AREA-SQUARE FEET 
TEMPERATURE-DEGREES FAHRENHEIT 
ENERGY- BTU PER HOUR 

IMPLICIT REAL*8( A-H.0-Z1 

THE NUMBER OF DATA POINTS, N, IS ENTERED. 

READ! 5, 5) N 
5 FORMAT (15) 

EPS=THE TOTAL WIND TUNNEL BOUNDARY CORRECTION FOR 
SURFACE BLOCKAGE AND WAKE BLOCKAGE. 

EPS=0 .0382L5 
G SUBC=3 2 . L 74 
R ANKN = 4 59 . 67 
PY =3 .14 15926535898 

STEFAN=THE STEFAN-BOLTZMANN CONSTANT. 

ST EFAN= 0.00000000 17 14 

D=D I AMET E R OF THE CYLINDER IN FEET. 

0=4.46875/12.0 

ALT 0=TH E TOTAL LENGTH OF NICHROME RIBBON IN THE MAIN 
CIRCUIT. 

A LT0=39 .0*PY*D+ 1.0/12.0 

R=RES ISTANCE OF THE MAIN HEAT GEN. CIRCUIT (OHMS). 
R=20.527 

CONDN I = THE THERMAL CONDUCTIVITY OF NICHROME. 

CONDN 1=7.5 l 

E MI SS = THE THERMAL EMISSIVITY OF NICHROME. 

EMISS=0 .7 

AC0ND=T HE AREA PERMITTING CIRCUMFERENTIAL CONDUCTION. 
AC ON 0=0.0031*0. 375/144 . 0 

THE VOLTAGE LEVEL, V, IS READ TO DETERMINE THE POWER. 
THE DYNAMIC PRESSURE, U,IS READ IN CM. OF WATER. 

• DT HET A=THE ARC LENGTH OF AN ELEMENT ASSOCIATED WITH 
ANY NODE/DATA POINT IN ANGULAR DEGREES. 

TU=THE TURBULENCE LEVEL OF THE WIND TUNNEL AT THIS 
FREE STREAM VELOCITY. 

RE AD ( 5, 10) V, U , DTHE TA , TU 
10 FORMAT (F 10.6, 10X , F 1 0. 3 , 10X , F 1 0. 3 , 1 OX ,F 1 0. 6 ) 

THE KINEMATIC VISCOSITY, CONDUCTIVITY 8 DENSITY OF AIR 
(VISC, COND, AND RO) ARE THE FILM TEMPERATURE VALUES. 
READ (5, 20 )VISC,COND,RO 
20 FORMAT( IX, 3( F 1 5 . 8 » 5X) ) 

TT CF=THE TOTAL WIND TUNNEL CORRECTION/CONVERSION 
FACTOR, AND IS DEFINED AS THE RATIO OF THE TUNNEL 
CORRECTION FACTOR, X, TO THE AREA CONTRACTION 
CONVERSION FACTOR, WHICH IS A FUNCTION OF A2/A1. 

TTCF= 1 . 0 

DELTAX = THE ARC LENGTH OF AN ELEMENT ASSOCIATED WITH A 
PARTICULAR NODE/DATA POINT IN FEET. 
DELTAX=D/2.0*DTHETA*PY/180. 0 
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ASURF=T HE SURFACE AREA OF AN ELEMENT ASSOCIATED WITH A 
PARTICULAR NODE/DATA POINT. 
ASURF=0.375*DELTAX/L2.0 
ASURFL =ASURF 

ALFA=CONDNI*ACOND/DELTAX 

QGEN=THE HEAT GENERATION RATE WITHIN A PARTICULAR 
ELEMENT, (BTU/HR) . 

QGEN=( ( V**2/R)/ALT0)*DELTAX*3600.0/L054. 35026449 
QGENL=QGEN 

THE VELOCITY IS EVALUATED IN FEET PER SECOND. 
VELU=DSQRT(U*2.0*GSUBC*TTCF/(R0*0.488256 196 ) ) 

THE VELOCITY IS CORRECTED FOR SURFACE, WAKE BLOCKAGES. 
VEL=D SQRT ( U*2 . 0*G SUBC* TTCF /( R0*0 . 488256 L 96 ) )*(L .0 + EPS) 
DYNHD=U *TTCF 

THE REYNOLDS NUMBER IS COMPUTED. 

REYNU=VELU*D*3600. O/VI SC 

THE REYNOLDS NUMBER IS COMPUTED USING CORRECTED VEL. 
REYN=VEL*D*36OO.0/VISC 

W RITE (6 , 22 ) 

22 FORMAT! • L' ,////////// ,45X,' FORCED CONVECTION HEAT TR', 
l'ANSFER FROM A UNI FORML Y • , / / , 5 I X , • HEAT ED CYLINDER IN', 
2* A CROSSFLOW OF AIR',///) 

WRITE (6 ,23) DYNHD, VELU,VEL, REYNU , REYN, D ELT AX , ASURF, 
IQGEN, TU 

23 FORMAT (44X, • FREE STREAM DYN . HEAD EQUALS ' , F6. 2, LX, 

I ' CM OF WATER',/ ,45X,'FREE STREAM VELOCITY EQUALS', F6. 
21, LX, ' FT .PER SEC . (UNCORRECTED) ' ,/ ,45X,'FREE STREAM ', 
3'VELOCITY EQUALS' ,F6.L, LX, 'FT. PER S EC .( CORRECTED )', / , 
446X , ' TH E REYNOLDS NUMBER E QUALS • , F8 . 0 , L X , ' ( UNCORR ECT • , 
5 ' ED ) ' , / , 46X , ' THE REYNOLDS NUMBER EQUAL S • , F 8. 0, LX , 

6* (CORRECTED) • ,/ ,59X,'DELTAX EQU ALS • , F8 . 5 , LX , ' F E ET ' , 

7/ ,60X, 'ASURF EQ UAL S • , F9. 6 , LX , • SQUARE FEET*,/ ,6LX, 
8'QGEN EQUALS' , F8 .5 , LX, ' BTU PER HOUR',/ ,63X,'TU EQUAL' 
9, 'S', F5. 2, LX, 'PER CENT',///) 

WRITE (6, 24) 

24 FORMAT (2LX,‘ THETA* ,2X, 'DELTA-P' , 2X , ' CP ' , 2X, 'DELTA-T' ,3 
LX, 'DELTA-T', 3X, 'DELTA-T', 4X,'QC0ND' ,4X,' QRAD' ,3X , ' HT * 
2 , 'TRANS' , 2X, ' NUSSELT ' , 2X, ' FROESSL ING' ,/ , 2 LX, • ( DEG ) • , 
33X, • (CM) • ,8X, ' (DEG F) ' , LX , ' CORRECT ION' , LX , • CORRECT ED* , 

4 LX, '( BTU/HR) • , LX, • ( BTU/HR ) • , 2X , 'COEFF' , 7 X , • NR ' , 8X , • NR' 

5 ,/) 

DO 50 I = L , N 

TDELM=T HE LOCAL TEMPERATURE DIFFERENCE BETWEEN THE 
NICHROME RIBBON AND THE AMBIENT TEMPERATURE IN MV. 
TDELT = T DELM CONVERTED TO DEGREES F. 

TAMM=T H E AMBIENT TEMPERATURE IN MILLIVOLTS. 

TAMB=TA MM CONVERTED TO DEGREES F. 

T FWM=T H E TEMPERATURE DIFFERENCE AT THE NEXT SUCCESSIVE 
NODE/DATA POINT IN MILLIVOLTS. 

TFWD = TF WM CONVERTED TO DEGREES F. 

T BACM=T HE TEMPERATURE DIFFERENCE AT THE PREVIOUS 
NODE/DATA POINT IN MILLIVOLTS. 

TBACK=TBACM CONVERTED TO DEGREES F. 

T CORM=T HE HEATED AIR AND VISCOUS HEATING CORRECTION TO 
DELTA T, IN MILLIVOLTS. 

T CORF=T CORM CONVERTED TO DEGREES F. 

THET A = T HE ANGLE MEASURED FROM THE FWD STAGNATION POINT 
READ( 5, 25) TDELM, TAMM, TFWM,TB AC M,TCORM, THETA 

25 FORMAT (2X,6(FL0.5, 2X) ) 

THE FOLLOWING EQUATIONS REPRESENT A COPP ER-CONSTANTAN 
THERMOCOUPLE CONVERSION TABLE. 

IF(TDELM.GT. (0.52L) )G0 TO 26 
TDELT=45.75L63398*(TDELM-0.2L5)+L0.0 
GO TO 28 

26 IF(TDELM.GT.( 0.854) )G0 TO 27 
TDELT=45.0L6077L7*( TDELM-0 . 543 ) + 25. 0 
GO TO 28 

27 TDELT= 44. L 6403785* ( TDE LM-O. 877)+40.0 

28 IF(TFWM.GT .(0.52L) )G0 TO 29 
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TF WD=45 .75163398* (TFWM-0 .2 15) +10 .0 
GO TO 3 1 

29 IF(TFWM.GT .(0.854) )G0 TO 30 

TF WD=4 5 . 01 60771 7* ( TFWM-0 .543) +2 5 .0 
GO TO 31 

30 TFWD=44. 16403785* (TFWM-0. 877) +40 .0 

31 IF(TBACM.GT. ( 0.521) )G0 TO 32 

T BACK= 45 .75163398*! TB ACM— 0 .215) + 10.0 
GO TO 34 

32 I F ( TBACM.GT. 0. 854) GO TO 33 

T BACK=4 5.0 16077 17 * (T B ACM-0 . 543 ) + 25 . 0 
GO TO 34 

33 TBACK=44. 1 6403785* ( TB ACM-0 . 877) +40.0 

34 TAMB=43 .2 19076* (TAMM-0 .765 ) +67.0 
TC0RF=46. 5116279*TC0RM 

C 

C T=T HE LOCAL TEMPERATURE AT A NODE/DATA POINT. 

T =TAMB+ TOE LT 
A SUR F=A SUR F 1 
QGEN=QGEN1 

C QCOND=THE NET HEAT FLOW INTO AN ELEMENT BY CONDUCTION. 

QC0ND=ALFA*(TFWD+TBACK-(2. 0*TDE LT) ) 

C QRAD=T H E HEAT FLOW RATE FROM AN ELEMENT BY RADIATION. 

QRAD= STEFAN* EMI SS*ASURF*( ( T+R AN KN) **4-( T AM B+R ANKN ) **4 ) 
C TDELCO=THE LOCAL TEMPERATURE DIFFERENCE CORRECTED FOR 

C THE NET EFFECT OF HEATED AIR AND VISCOUS hEATING. 

TDELCO = TDE LT-TCORF 

C HCONV = T HE LOCAL FORCED CONVECTION HEAT TRANSFER COEFF. 

HCONV=( QCOND+CGEN-QP.AC) / ( ASURF*T DELCO ) 

C GNU=THE LOCAL NUSSELT NUMBER. 

GNU=HCONV*D/COND 

C THE FROESSLING NUMBER IS CALCULATED. 

F=GNU/DSQRT(REYN ) 

C THE LOCAL DELTA P IS READ IN INCHES OF WATER. 

R EAD( 5*43) DELTAP 
43 FORMAT ( F10 .3) 

C THE PRESSURE COEFFICIENT IS CALCULATED. 

DELTAP = 2. 54*DE LTAP 
CP=DELTAP/DYNHD 
C 

WR IT E( 6, 49) THETA, DELTA P, CP, TDEL T , TCORF , TD ELCO , QCOND , 
1QRAD,HC0NV ,GNU , F 

49 FORMAT! 19X ,F7 . 1 ,2X , F6 .2 ,2X , F5 .2 , IX , F6 .2 , 2X , F7 .2, 3X , F7 . 
12,3X,F 8.5, IX, F8.5, 2X, F6 . 2 , 3X ,F 7 . 2 , 2X ,F 8 . 4) 

50 CONTINUE 
STOP 

END 
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APPENDIX C 



Wind Tunnel Calibration for the Combined Effect 
of Ambient Air Temperature Rise and Viscous Heating 

of the Cylinder 

In general the free stream air temperature steadily 
increased during the operation of the wind tunnel. This 
temperature rise was attributed to the viscous heating in 
the boundary layer of the wind tunnel walls and the energy 
added by the rotating impeller. At low speeds (Re<260 ,000 )' 
this temperature rise was very gradual and did not affect 
the' accurate measurement of AT. The measured AT represented 
the temperature rise solely due to the applied power. Al- 
though the surface of the cylinder became warmer due to the 
heated air this extraneous heating of both. the cylinder and 
the air were at the same rate, and AT remained a true 
measurement. During high speed runs, however, the rate 
of temperature rise of the tunnel air was rapid and ex- 
ceeded the corresponding cylinder surface temperature rise. 
This effect caused the measured AT to be smaller than it 
should have been. 

At the same time, and also for the high speed runs, the 
viscous heating of the cylinder surface provided a counter 
effect, tending to increase the measured AT. In all of 
the high speed runs this latter effect was dominated by the 
ambient air temperature rise. 
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Rather than attempting to uncouple these two effects, 
a scheme was devised to calibrate their combined effect 
out of the data for the high speed runs. The Nichrome rib- 
bon was not energized and the cylinder was heated only by 
the extraneous effects discussed in the preceding para- 
graphs. After several hours steady state conditions were 
determined using the same technique as described in Section 
V. EXPERIMENTAL PROCEDURE for an ordinary run. The temper- 
ature difference ( T 0 -T a j_ r ) was then measured at various 
angular locations around the cylinder, and treated as a 
correction factor for the data collected during an actual 
run with applied power at that Reynolds number. 

Calibration runs were performed at Reynolds numbers of 
309,000, 406,000, and 495,000, and the AT correction factors 
are provided in Section VI-A TABULAR RESULTS. In the data 
reduction computer program (Appendix B) the AT correction 
factors are subtracted from the measured AT just prior to 
the calculation of the local heat transfer coefficient. 

The radiation and conduction losses are based upon the un- 
corrected, measured AT. The AT correction factors are most 
often negative because the heated air effect was dominant. 
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APPENDIX D 



Wind Tunnel Turbulence Level Measurement 

The turbulence level of the Aerolab wind tunnel test 
section was measured with a Thermo-Systems hot wire anemo- 
meter at eleven different free stream velocities ranging 
from 10 to 110 mph. 

Initially, a Thermo-Systems calibrator model 1125 was 
used to calibrate the probe. The cold wire (room tempera- 
ture) resistance was 7*77 ohms. With zero air flow through 
the calibrator the initial probe voltage was 1.667 volts. 
The probe voltage was then recorded for various controlled 
values of air flow velocity through the calibrator. The 
results of this calibration are presented in Fig. 18, as 
a plot of the probe voltage squared, V 2 , versus the square 
root of the air flow velocity through the calibrator, /U. 
The intercept of the linear calibration curve which result- 
ed, was the square of the actual value of initial probe 
voltage, V Q , to be used in the turbulence level 
calculations . 

The probe was then mounted in the center of the test 
section on a vertical rod whose height could be adjusted 
in one inch increments. At each air speed for which turbu- 
lence level measurements were taken, the rod was used to 
traverse the vertical height of the test section. The 
probe voltage, V, and the rms voltage, e, were recorded at 




/u 
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five inch intervals along the vertical. These readings 
were observed to be constant for any given air speed except 
near the floor and ceiling where wall boundary layer ef- 
fects caused some fluctuation. 

The turbulence level, Tu, was calculated at the various 
air speeds by making use of the relation 

Tu vjz = ijlv 

u v 2 -v 2 

o 

where U represents the mean free stream velocity, and is 

the mean fluctuation in the free stream velocity. The 
results are shown in Pig. 19 . 
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-0 10 20 30 40 50 60 70 80 90 100 110 

Free Stream Velocity, U, (mph) 

Figure 19. Turbulence Intensity Level of the Wind Tunnel Determined at Various Air 
Speeds 



APPENDIX E 



Endurance Test for the Acrylic Tube 

An eighteen in section of the 4.5 in o.d. cast acrylic 
tubing was wrapped with six turns of the Nichrome ribbon. 
The ribbon was instrumented with one copper-constantan 
thermocouple. The surface temperature was then measured 
using a Leeds and Northrup millivolt potentiometer and an 
ice bath for a reference. A Lambda 60 volt power supply 
was used to provide a 2. 0-2. 5 amp current through the 
ribbon. The tube was thereby heated, uniformly and slowly, 
to temperature levels of 100°P, 110°P, 120°P and 150°P, 
and cooled slowly to room temperature after each heating 
cycle. The heating and cooling rates were only about one 
to two degrees P per min., to minimize the relative differ- 
ences between the acrylic and the Nichrome during the ex- 
pansion and contraction. In all of the heating cycles 
except the final one the cylinder did not physically de- 
form,' as there was not a noticeable expansion of the 
acrylic. In the final test, between 130-135°P the acrylic 
began to visibly expand causing the ribbon to tighten. 

This immediately caused a hoop stress on the surface which 
deformed the cylinder by constraining its expansion in the 
area where the ribbon was wrapped. The cylinder expanded 
further in the unwrapped areas. When the cylinder cooled 
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it contracted uniformly leaving a smaller, deformed diameter 
in the region of the ribbon. 

Prom this test it was concluded that to avoid the risk 
of deformation the actual test cylinder should not be 
heated above 120°F. 

In a second test the surface temperature was heated to 
l80°F, severely deforming the cylinder. However, the 
R.T.V. sealant was unaffected and held solidly. 



/ 
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APPENDIX F 



Uncertainty Analysis 

As in all experimental studies, there existed some 
degree of uncertainty with regard to the accuracy of vari- 
ous results, which stemmed from the data collected. To 
evaluate the levels of uncertainty involved, the method of 
Kline and McClintock described in Ref. 44 was employed. 

The measured variables which were the origins of the 
uncertainty were APj_ 2 , AP ^ n d er , AT, V, and R (see no- 
menclature for the definitions of these terms). The fluid 
properties of air were extracted from table A-7 of Ref. 45, 
and their uncertainty is neglected in the analysis. 

The magnitude of the uncertainty associated with the 
pressure measurements depended upon the flow velocity of 
the run, and the angular location on the cylinder. Simi- 
larly, the magnitude of the uncertainty associated with the 
temperature measurements depended upon the angular location. 

In this appendix the calculations of uncertainty are 
shown in detail for a representative Reynolds number of 
260,600, and at the forward stagnation point of the cylin- 
der. Uncertainties for other flow velocities and locations 
are discussed. 

The Pressure Coefficient 

The pressure coefficient was calculated from equation 

28 , 
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c 



p 



AP 

cylinder 

APl-2 



Applying the technique referenced above, one has 



1 

a* 

o 

3 


,(0 AP 

cylinder 


2 

+ 


[ W A Pl _ 2 j 


o 

hd 

i 

S 


k AP cylinder 


X 


A?i -2 



where co is the uncertainty interval, or precision range. 
At the forward stagnation point, 

AP cylinder = 7 * 5 - °‘ 1 cm water (10 to 1) 

APj_ 2 = 7*7 ± 0.1 cm water (10 to 1) 

Substituting, one has 



0 ), 



0.10 


2 

+ 


'cT 

1 — 1 
0 

V 


17.52J 


T 


17.71J 



« 0.02 



Therefore, the uncertainty in the pressure coefficient was 
about two per cent. In the wake, the uncertainty was about 
four per cent because the unsteady motion of the air les- 
sened the readability of the manometer. 

At the higher Reynolds numbers the uncertainty in the 
pressure coefficient was about three per cent at the for- 
ward stagnation point and about six per cent in the wake. 
The higher speed flows produced larger fluctuations on the 
manometer. 

At the lower Reynolds numbers the measured pressures 
were small and the uncertainty intervals were a significant 
fraction of these pressures. As a result the uncertainty 
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In the pressure coefficient was large, reaching as high as 
49% for the lowest run. (Re = 57,000) 

The Heat Transfer Coefficient 

The heat transfer coefficient was calculated from the 
equation, 

V 2 

h - R 

n 6 A AT 
s 

where A g is a constant (see nomenclature for definition of 
terms). Continuing, 



II 

■Cl 
3 1 
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> 
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C\J 

V. 
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/ 

W R 
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“at 


h / 


v J 




L R J 




AT 



For the experimental run at a Reynolds number of 260,600, 

V = 59.1 ± 0.3 volts (10 to 1) 

R = 20.5 ± 0.2 ohms (10 to 1) 

AT = 20.0 ± 0.5°F (10 to 1), at the forward stagnation 
point . 

Substituting, one has. 



“h / 


[ 2 o o- 3 ] 


2 

4. 


° • 2) 2 + 


f 0.5] 


h / 


r 59. ij 


1 


[ 20 . 5j 


120. Oj 



The uncertainty in the heat transfer coefficient was about 
three per cent at the forward stagnation point for all the 
experimental runs. 

At the transition point where the boundary layer be- 
comes turbulent, or at the separation point in laminar 
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flow, the temperature difference had an uncertainty Inter- 
val of about five°F, and the uncertainty in the heat trans- 
fer coefficient was on the order of 25$. 

The uncertainty of the heat transfer coefficient in 
the wake was about six per cent. 

The Nusselt Number 

The Nusselt number was calculated from equation 3, 




Continuing , 




and 

D = ± .05 in (10 to 1). 

The uncertainty of the measurement of the diameter was 
negligible. Therefore, the uncertainty in the Nusselt 
number was the same as the uncertainty in the heat transfer 
coefficient. 

The Reynolds Number 

The Reynolds number was calculated from equation 1, 

Ro - UD 
Re " — ’ 

f 

and its uncertainty is the same as that of the free stream 
velocity, U, which was calculated from equation 29, 
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Continuing , 



ll)„ (Or, 

Re _ U 
Re “ U 



0.5 



(0 



AP. 

AP 



1-2 

1-2 



0.5 



O.l' 
J .71 



0.007. 



Therefore, the uncertainty in the Reynolds number is approx- 
imately half the uncertainty in the dynamic pressure, and 
equals about 0.7 %• At the high Reynolds numbers this un- 
certainty doubles. 



The Froessling Number 

The Froessling number was calculated from equation 4, 

F = M- . 

/Re" 

Continuing , 
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[2ReJ 



Substituting from previous results, one has 

-pr = / ( 0 . 03 ) 2 + (0 . 0035 ) 2 * 0.03 • 

The uncertainty of the Froessling number was generally 
about the same as that of the Nusselt number. At the for- 
ward stagnation point the uncertainty of the Froessling 
number was about 3 -^% in most cases. 
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